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Abstract 



The analysing power measurements for the pp pprj reaction studied in this 
dissertation are used in the determination of the reaction mechanism of the 7] 
meson production in nucleon-nucleon colhsions. 

Measurements have been performed in the close-to-threshold energy region at 
beam momenta of Pbeam — 2.010 and 2.085 GeV/c, corresponding to the excess 
energies of Q = 10 and 36 MeV, respectively. The experiments were realised by 
means of a cooler synchrotron and storage ring COSY along with a cluster jet 
target. For registration of the reaction products the COSY-11 facility has been 
used. The identification of the rj meson has been performed with the missing mass 
method. 

The results for the angular dependence of the analysing power combined with 
the hitherto determined isospin dependence of the total cross section for the rj 
meson production in the nucleon-nucleon collisions, reveal a statistically signifi- 
cant indication that the excitation of the nucleon to the Sn resonance, the process 
which intermediates the production of the r) meson, is predominantly due to the 
exchange of a vr meson between the colliding nucleons. 

The determined values of the analysing power at both excess energies are 
consistent with zero implying that the f] meson is produced predominantly in the 
s-wave at both excess energies. 

Streszczenie 

Zaprezentowane w tej pracy pomiary zdolnosci analizuj§cej dla reakcji 
~^ PPV majg, na celu wyznaczenie mechanizmu produkcji mezonu r) w zderzeni- 
ach nukleonow. 

Pomiary wykonano w przyprogowym obszarze energii. W ekspery- 
mentach wykorzystano wi§,zk§ protonow spolaryzowanych poprzecznie o 
p§dach Pbeam — 2.010 i 2.085 GeV/c, odpowiadajcicych energiom wzbudzenia 
Q = 10 oraz 36 MeV dla reakcji ;pp — > ppr). Do badan wykorzystano synchrotron 
COSY, tarcz§ klastrow§,, oraz system detekcyjny COSY-11. W doswiadczeniach 
rejestrowano produkty reakcji obdarzone ladunkiem elektrycznym, natomiast 
identyfikacji mezonu rj dokonano za pomoc§, metody masy brakuj§cej. 

Uzyskane doswiadczalnie rozklady k§,towe zdolnosci analizuj§,cej, w pol§,czeniu 
ze zmierzong, uprzednio zaleznoscig, izospinowg, produkcji mezonu 77 w zderzeni- 
ach nukleon-nukleon, pozwalaj§, stwierdzic, iz wzbudzenie nukleonu do rezonansu 
Sii(1535) - proces posredniczacy w produkcji mezonu r] - nastepuje w glownej 
mierze poprzez wymian§ mezonu tt pomi§dzy zderzaj§-cymi si§ nukleonami. 

Wartosci zdolnosci analizujcicej dla obu energii wzbudzenia sg, rowne zeru w 
granicach bl§d6w statystycznych, co wskazuje, iz mezon 77 jest w glownej mierze 
produkowany w fali s. 
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Chapter 1 
Introduction 



The main goal of this dissertation is the closer insight into the production mech- 
anism of the rj meson in the interaction of nucleons. Despite the fact that the 
discovery of this meson - a member of the pseudoscalar meson nonet (see Ap- 
pendix |A]) - took place almost half a century ago [1] its production dynamics has 
remained an open question for a long time. 

From the precise and extensive measurements of the total cross section for the t] 
meson production in the pp — ^ ppr] reaction in the close-to-threshold region [2-10] 
it was concluded [11] that this process proceeds through the excitation of one of 
the protons to the Sii(1535) state which subsequently deexcites via the emission 
of the 7] meson and a proton. The crucial observations here were the large value 
of the absolute cross section and isotropic distributions of the angle of the t] 
meson emission in the reactions centre-of-mass system [8, 12]. However, there 
are plenty of possible scenarios of the excitation of Sii(1535). In fact, any of 
the TT, 1], uj or/and p mesons may contribute to the resonance creation. Taking 
solely into account the total cross section one cannot deduce which one out of this 
rich spectrum of mesons plays the most important role in the excitation of the 
intermediate resonance. 

Further investigations of the rj meson production process, namely the 
determination of the isospin dependence of the total cross section by the 
WASA/PROMICE collaboration [13], put some restrictions to the above- 
mentioned possibilities of the resonance excitation process. The ratio 

= aijm — i> pnri)/a{pp — * ppr]) was determined to be about 6.5 in the excess 
energy range between 16 and 109 MeV, which revealed strong isospin dependence 
of the production process. The production of the 77 meson with the total isospin 
1 = exceeds the production with the isospin I = 1 by a factor of 12, suggest- 
ing [13] that the isovector meson exchange - the vr or/and p meson exchange - is 
the dominant process in the excitation of the Sn resonance. However, the relative 
contributions of the pseudoscalar tt meson and vector p meson still remained to 
be determined. 

Here the measurements of the polarisation observables can assist, because the 
predictions of the one boson exchange models [11, 14] with respect to the most 
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basic polarisation observable as the analysing power are sensitive to the type of 
meson being exchanged between the protons in order to excite one of them into 
the resonant state. Measurements of the beam analysing power for the frp — > pprj 
reaction performed by the DISTO collaboration [15] as well as the tentative ex- 
periment of the COSY-11 group [16, 17] did not bring the univocal solution of this 
problem. The interpretation of results of the DISTO measurements, performed in 
the far-from-threshold region at the excess energies of Q = 324, 412, and 554 MeV 
suffered from the lack of a theoretical prediction for the analysing power. This 
is due to the fact that far from the reaction threshold the higher partial waves 
are involved in the reaction process, and so the theoretical description becomes 
more complicated. COSY-11 results at the excess energy of Q = 40 MeV, due to 
the insufficient statistics, could not have been used in order to judge between the 
predictions of pseudoscalar and vector meson exchange models. 

This gap is filled by the results of experiments described in this dissertation. 
For the first time ever the COSY-11 group was able to perform the measurements 
of the analysing power for the pprj reaction in the close-to-threshold region 
that has brought the answer to the open question of the rj meson production 
mechanism in hadronic collisions. 

It is worth mentioning that the precise determination of this mechanism - a 
process that proceeds through the strong interaction - delivers information about 
the nature of the strong forces. Thanks to the results presented in this dissertation 
the theoretical models might be revisited with new input parameters applied to 
these models: the coupling constants in the description of the production process 
of the 1] meson, the initial and final state interactions and also dimensions of the 
reaction region. 

Experiments presented in this work have been performed by the COSY-11 
collaboration by means of the COSY-11 facihty [18-20] at the COoler SYnchrotron 
and storage ring COSY [21,22] in the Research Center Jiilich in Germany. The 

analysing powers have been measured during two runs at difi^erent beam momenta: 
Pbearn = 2.010 GcV/c (May 2003) and 2.085 GeV/c (September 2002), which for 
the pp ppr) reaction correspond to the excess energies of Q = 10 and 36 MeV, 
respectively. 

This thesis consists of eight chapters. In the following one the basic definitions 
are given, which are important in the understanding of the general methods that 
are applied to the experiments with polarised beams. In particular the definitions 
of the spin of a quantum particle, the polarisation vector, and the analysing power 
are given. We introduce the reference frames used in the polarisation experiments, 
in particular we show how to construct the so called Madison frame following 
the Madison convention. This is relevant in order to determine the sign of the 
analysing power. At the end of this part of the thesis we derive the formula which 
will be used to determine the analysing power. 



3 



In Chapter [3] an overview of the theoretical models of the t] meson production 
in hadronic collisions is sketched. Superimposed are the phenomenological models 
working within the one boson exchange frame. A special emphasis is put on the 
pseudoscalar meson exchange model of Nakayama et al. [11] and the vector meson 
exchange model of Faldt and Wilkin [14] as these are the only models giving 
the predictions of the analysing power for the pp ppt] process in the close- 
to-threshold region. The summary of the so far measured observables for the 
pp — > ppt] reaction and the conclusions from these measurements are given in this 
chapter. 

The general description of the cooler synchrotron COSY along with the method 
of the formation and acceleration of the polarised proton beam are elucidated 
in Chapter |H Also the experimental facility COSY-11 is briefly described. The 
detector setup and the method of analysis has been previously presented in various 
dissertations [23-25] so here we will only focus on the description of the detectors 
which are interesting from the point of view of the analysis of the pp — ^ ppi] 
reaction. 

In Chapter O the calibration of the detectors needed to register the pp — *■ 
ppi] process is described. Also the methodology of the analysis is given i.e. the 
missing mass technique, identification of the beam spin mode from the number 
of elastically scattered events, the methods of determination of the production 
yields, relative luminosity, and the beam polarisation. At the end of this chapter 
the final results of both analyses are presented. 

Chapter [6] contains the confrontation of the theoretical models with the exper- 
imentally determined values of the analysing power. Using statistical inference we 
conclude about the production mechanism of the rj meson in hadronic collisions. 

The possibilities of extending the experiments presented in this work are dis- 
cussed in Chapter [71 

In Chapter [8] we summarize the thesis. Conclusions from the analyses pre- 
sented in this dissertation can be found in this part. 

Three appendices have been added at the end of this dissertation in order to 
elucidate in more detail some issues discussed in the text. We shall explain the 
differences between the pseudoscalar and vector mesons (Appendix , prove a 
property of the analysing power (Appendix [B]) postulated and used in the analysis 
in Chapter [61 We will also demonstrate that the elastic scattering in the polarisa- 
tion plane does not depend on the degree of polarisation (Appendix [C|) , a property 
which originates from the parity invariance rule and is used in Section 15.31 for the 
calculation of the relative luminosity. 



Chapter 2 
Definitions 



2.1 Spin and the beam polarisation 



The formal quantum mechanical description of spin 1/2 particle is given. The polarisa- 
tion vector and polarisation plane are defined. 

The spin of a particle is an internal degree of freedom, which is governed by the 
same equations of motion as the angular momentum. Therefore, we can assign to 
observable "spin" a vectorial operator S = [Sx, Sy, Sz], where Si (i=x,y,z) are the 
hermitian operators. Subsequently we will only consider the spin- 1/2 particles. 
In this case the spin operators can be expressed in terms of the Pauli matrices af. 

Si = ^cji, (2.1) 

where 

? J) (i (° o') • 

With respect to the arbitrary quantization axis Oy, a spin 1/2 particle can be 
found in one of the two eigenstates of Sy. These eigenstates may be represented 
by the eigenfunctions 

The x± are the eigenfunctions of Sy with the eigenvalues ±h/2. In what follows 
we will refer to the x+ X- states as to the "spin up" and "spin down" states, 
respectively. 

Let us denote the fractions of protons in the pure spin up and spin down states 
by n+ and n_, respectively. Define a polarisation Py along the Oy quantization 
axis as the asymmetry of the populations n+ and n_ : 

Py = ^ . 2.4 
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Analogously we could define the polarisations and with respect to the Ox 
and Oz axes, but we will assume that these two cancel out in the experiment 
described in this thesis. From these properties we can construct a polarisation 
vector P defined as follows: 

P = [P^,Py,P,] = [0,P,0]. (2.5) 

The modulus of P shall be called the degree of polarisation. A single spin 1/2 
particle in a pure spin state x+ or X- is fully polarised, and for such a particle 
I P I = 1 , which follows from definition 12.41 and the fact that single fermions can 
be found only in one out of two states: x+ o^' X-- lu reality we deal with the 
systems containing a large number of particles, and for such systems the degree 
of polarisation is usually smaller than 1. 

By the polarisation plane we will refer to a plane spanned by the beam mo- 
mentum vector Pbeam = [0, 0,Pz] and a polarisation vector P = [0, P, 0]. 

2.2 Reference frames 

The reference frames along with the angles of production are defined. 

We will describe the polarisation observables within the Cartesian coordinates, 
however in the beginning we have to introduce the certain reference frames that 
are convenient to this description. 





' M 

y 




Pbeam 




M 





Figure 2.1: Definition of the reference system according to the Madison convention 
for the description of the polarisation observables. 

According to the Madison convention [26] we construct the reference system 
(x'^, y'^, z^) in the following manner: let a unit vector of the 02;*^-axis be 



2.2 Reference frames 
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parallel to the beam momentum direction given by the vector Pbeam, and the 
0|/*^-axis unit vector y'^ along the direction of the Pbeam x Pout vector, where 
Pout is the momentum vector of the i] meson. Both pbeam and pout vectors 
are defined in the centre-of-mass (CM) frame. The x'^ - a unit vector along 
the Ox*^-axis is defined such that the basis (x'^,y^,z^) forms an orthonormal, 
right-handed set. Such defined reference frame is sketched in Fig. 12. 1[ We will 
refer to this frame as to the Madison frame. 

The angle 9^ will be called the polar angle of rj production. This is an angle 
between pout and z^: 



cos 



Pout ■ Z 



M 



(2.6) 



out 



Now, we introduce the fixed in the centre-of-mass an accelerator coordinate 



frame (x 
momentum 



y ) ^ ) 



Define z***^*^ - a unit vector that is parallel to the beam 
, a unit vector y*^*^*^ perpendicular to the accelerator plane, pointing 
up, and x*^*^*^ which completes the right-handed basis. The accelerator frame is 
shown in Figure [221a. 




Pbeam 



= 71, right 



Figure 2.2: (a) Accelerator reference frame. Definition of the scattering to the left 
(b) and to the right (c) with respect to the polarisation plane. 

In the case of the accelerator frame the polar angle of r/ - - is the same as 
in the case of Madison frame, as the Oz'^^ and Oz"''^'^ axes are defined in the same 
manner. 

By the azimuthal angle of the rj meson emission - 0,^ - we will understand 
the angle between the x**^^ vector and the projection of Pout onto the x"^^-y"'^'^ 
plane. Both angles are presented in Figure I2.2[ a. If the momentum of the rj 
meson is lying in the horizontal plane, i.e. if 0^ = or 0.^ = vr, then the Oy^^ and 
Qyacc g^jjgg (and also Ox*^ and Ox"''^^ axes) may be either parallel or antiparallel, 
depending on whether the rj meson is emitted to the left (Fig. I2.2[ b) or to the 
right side (Fig. I2.2[ c) in the accelerator system. 

It is important to distinguish between the Madison and the accelerator frame. 
Whereas the accelerator frame is fixed in the centre-of-mass frame, the Madi- 
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son coordinate system may vary in space from event to event. According to tlie 
Madison convention all the physical observables, like for example the beam polar- 
isation vector, should be considered in the Madison frame in order to avoid the 
ambiguities in the determination of the sign of the analysing power. 



2.3 Analysing power 



A five-dimensional phase space for the description of the A^A^ NNM reactions is 
defined. The most general formula for the cross section for experiments with polarised 
beam and target is given, which contains the definition of the analysing power. 

For a nuclear reaction with a given initial channel and a three body final state 
twelve parameters (four-momenta of three particles) have to be known to fully 
describe the exit channel. Some of these twelve parameters are bound with the 
relativistic formula: 

^' = ^' + (Pi)', ^ = 1,2,3; (2.7) 

where Ei, rrii, pi denote the total energy, mass, and the momentum vector of the 
i-th particle. Equations 12.71 reduce the number of variables to nine, and these 
variables are dependent on the initial state parameters on the basis of the four- 
momentum conservation (four additional equations). Therefore, there are eventu- 
ally five linearly independent variables that need to be measured for each reaction 
in order to fully describe its kinematics. Here, we follow the references [8, 27] in 
the choice of these kinematical variables. We will use the invariant masses of the 
proton-proton system rupp and proton-r/ system rUpri, the polar 6n and azimuthal 
(prj angles of the rj momentum in the centre-of-mass (CM) frame, as well as the an- 
gle ip - describing the rotation around the direction established by the momentum 
of the rj meson [27]. These variables are orthogonal and form a basis in the five- 
dimensional phase space. Let this basis be denoted by C = {^pp, ^prj, 4>r]^ ^v^ 
Let us denote by 

'^^'^^ dnippdmpr^dipridOridip ^ ^ 

the differential cross section for the reaction in absence of the polarisation of beam 
and target, and since the internucleon strong interaction is spin-dependent, let us 
also define 

a(C, P, Q) ^ ^MCP.Q) .2.9) 

' ' dmppdnipridcpndOridtp 

as the differential cross section for the reaction induced by a polarised beam on the 
polarised target, where P and Q are the polarisation vectors of beam and target, 
respectively. The beam polarisation vector P and the target polarisation 
vector Q should be referred to in the Madison frame defined in Section [2l2l 



2.3 Analysing power 
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In the most general case, when both P and Q are the non-zero vectors, the formula 
for cr reads: 

3 3 3 

a{C, P, Q) = (To(C)(l + P^MO + E ^'^'^(0 + E ^^^(P' (2.10) 

i=l i=l ')i=l 

where Aj and A^ are the beam and target analysing powers, respectively, and 
Cij{() are the spin correlation coefficients. Pi, Qi and Cij are the expansion coef- 
ficients, which depend on the polarisation of beam and target. In our case, where 
the target is unpolarised, the Qi and Cij coefficients disappear and we only deal 
with the Ai vector. Hence Equation 12.101 simplifies to: 

3 

a{C,P)=cro{0{l + Y,PMO)- (2.11) 

i=l 

For in Equation 12.111 the scalar product of P and A is present, in the case 
of a transversally polarised beam we can rewrite this equation in the following 
manner: 

a(C, P) = ao(C)(l + PAyirUpp, rupr,, 6^,, i/j) cos(0^)), (2.12) 

where P is the component of the beam polarisation normal to the incident beam, 
such that the ^/-component of polarisation vector expressed in the Madison frame 
is given by Py = Pcos0^. 

At this point we would like to mention that the acceptance of the COSY-11 
facility allows to register only events scattered near the horizontal plane. In the 
analysis the azimuthal angle 0,, was restricted to values of cos (j)ri ranging between 
0.87 and 1. This is depicted in Figure 12.31 where we show the shape of the 
acceptance of the COSY-11 system for the registration of the events with the 
Tj meson production as a function of the azimuthal angle 0^. One can see that 
the acceptance for events with 0.^ G [30°, 330°] is relatively very small. Moreover, 
the acceptance is much larger for the cases where the r] meson is emitted to the 
left (0.^ ~ 0°) than when it is emitted to the right side of the polarisation plane 
{(pn ~ 180°). Therefore we decided to consider in our analyses only the events 
with the 7] meson production to the left side with respect to the polarisation plane. 
For such a case cos 0,^ ~ 1, and considering that the beam was vertically polarised 
along the Oy"'^'^ axis (polarisation vector given by Equation 12. 5p we can rewrite 
Formula EH 

e,,^)). (2.13) 
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I 

50 100 150 200 250 300 350 

Figure 2.3: Distribution of the pp — > ppr] events wliich can be registered witli the 
COSY-11 setup as a function of the azimuthal angle of the rj meson production in the 
centre-of-mass system - 

It is worth mentioning that in the case of the pp — >■ ppr] reaction A;j.{() and 
Ay{() are equivalent, because the rotation of the whole system by tt/2 around the 
Oz axis does not lead to a change of the scattering pattern. Therefore, in case of 
spin 1/2 particles, there are only two independent polarisation observables: Ay{() 
and Az{C) - the transversal and longitudinal analysing power, the latter being of 
no significance in the experiments considered here. 

From Equation 12.131 it follows that the vector analysing power Ay{() may be 
understood as a measure of the relative deviation between the differential cross 
section for the experiments with and without polarised beam (in the absence of 
the target polarisation), normalized to the beam polarisation: 

MO-^'-^^^^f^. (2.14) 

2.4 Practical formula for the analysing power 

The practical recipe for the calculation of the analysing power is derived. 

In this section we consider the meson production where the beam consists 
of vertically polarised protons and the production takes place on an unpolarised 
proton target. 



2.4 Practical formula for the analysing power 
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Naturally, the Madison frame defined in Section is not fixed in space, as is 
schematically shown in Figure I2.4[ In the left panel there is depicted a situation 
when the rj meson is produced to the left side with respect to the polarisation 
plane, whereas the right panel shows the production of the t] meson to the right 
side with respect to the polarisation plane, defined in Section 12. 1[ While the 
Madison frame changes depending on pit the accelerator frame remains unalterred 
in space. One should also notice that the production of the r] meson to 
the left (right) with the polarisation vector of the beam pointing up in 
the accelerator frame is physically equivalent to the production of the 
T] meson to the right (left) with the beam polarisation vector pointing 
down. In both cases the physics, meaning the spin-dependent nuclear interaction, 
is the same and the two abovementioned cases are equivalent, from which we will 
make use in due course. 




O 



a) 





M acc 
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Figure 2.4: Schematic picture of the r] meson production to the left (a) and right side 
(b) with respect to the polarisation plane. The polarisation vector P is pointing up (out 
of the paper), which is denoted by Q. Black squares are the virtual rj "detectors". 



In Figure 12.41 we presented schematically the production of the rj meson in the 
accelerator frame's plane x"''^'^-z"''^'^ . The polarisation vector of the beam of protons 
is pointing up, i.e. along the y"'^'^ axis. Denoting by N(9n, 0,^)0 the number of the 
rj mesons emitted into the solid angle around 6'^ and 0^ and by Ni and N2 the 
numbers of the rj mesons reaching the virtual 1] "detectors" 1 and 2 we can write: 

Nl = N^{e„0) = Nl = cxo{er,){l 

and 

Nl=N^e„'K) = Nl = ao{e^){l 



We remind, that the angles 0,, and </>,, are referred to the accelerator frame. 



PU,(e,))E(e„0) j LUt, (2.15) 



P^Ay{e^,))E{er,,'n) / LUt, (2.16) 
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where the symbol "|" stands for the spin up orientation, J L^dt is the luminosity 
during the measurement with spin up, integrated over the whole time of the 
measurement, -^(6*^, (pr)) is the efficiency of the detection setup for "registering" 
the Tj mesons produced into the solid angle around Or^ and 0,,, aoiOr,) is the cross 
section for the rj meson production with absence of the beam polarisation, 
is the beam polarisation during the spin up cycles and AyiO^^) is the analysing 
power for the rj production averaged over rUpp, rUpri and ip. The "+" sign denotes 
the production with the polarisation vector along the Oy'^-axis (according to 
the Madison convention Py is positive in this case), and the "— " sign is for the 
production with the polarisation vector antiparallel to the Oy^^-axis (negative 
Py)- 

Now, if we flipped the polarisation vector i.e. if we made it pointing down 
(into the paper in Figure [23]) we would get: 

Ni = N^{er„o) = Ni = aoie,)ii-p^Ay{e,))E{e„o) j lUt, (2.17) 

= iVi(^„7r) =Ni = ao(e,)(l + pi/l,(^,))E(^„7r) j L^dt, (2.18) 

where the quantities with arrows pointing down refer to the respective quantities 
from Equations 12.151 and I2.16[ but for the cycles with spin down mode. 
From Equations 12.151 and 12.161 we obtain: 

The same is valid for Equations 12.171 and I2.18[ but we have to change t into | in 
Equation I2.19[ Note, that Equation 12.191 is independent of the integrated lumi- 
nosity. Equation 12.191 shows that the analysing power Ay{9rj) may be calculated 
from the production with only one spin orientation (spin up or spin down), if 
both scatterings - into left and right side in the accelerator plane - can be mea- 
sured. To some extend it may be realised by the COSY-11 detector. However, 
for the COSY-11 detector setup E{9ri,0) ^ E{9rj,n) (see Figure [231 where be- 
tween (pri G [330°, 30°] there are much more events than for 0^ G [150°, 210°]), and 
therefore in practice we have to perform the measurement with two beam spin 
orientations. 

Assuming that the degree of polarisation during the cycles with spin down is 
equal to the degree of polarisation for cycles with spin up, which is the case within 
the 2% of accurac}!! as has been shown in the previous measurements by means of 
the EDDA polarimeter [28], we may introduce the average degree of polarisation 
P: 

pT _|_ pi 

P ^ P^ ^ P^ ^ . (2.20) 



^See also Section [5X2l 



2.4 Practical formula for the analysing power 
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Assuming that the degrees of polarisation for spin up and down are equal, the 

T I 

asymmetry between and may be used as a measure of the spin-dependent 
r) meson production. Dividing Equation 12.151 by Equation 12.171 and replacing 
and by P yields to: 



NliOr^) _ J L^dt l + PA 



y\"vJ 



Introducing the relative luminosity defined as: 

/ L^dt 



-'rel 



and solving Equation 12.211 for Ay we obtain: 



(2.21) 



(2.22) 



Ay{9„) = re^^^ 2.23 

'''' p Nlie,) + u^Mie,) 

a formula which is independent of the efficiency of the detection setup. In what 
follows, Formula 12.231 will be used for the calculation of the analysing power. 



Chapter 3 
Theory 



3.1 Theoretical models of the 77 meson produc- 
tion in nucleon-nucleon collisions 

Some phenomenological models of the r] meson production in the pp ppr] reaction are 
shortly reviewed. Conclusions from the analysis of the available physical observables 
along with the predictions for the analysing power are pointed out. 

Despite the fact that the discovery of the rj meson took place over forty years 
ago [1], its production mechanism still remains an open question [29,30]. Af- 
ter early measurements of the total cross sections for the pp ppt] reaction in 
bubble chamber experiments [31-38], only recently there appeared high-statistics, 
close-to-threshold data from storage rings, giving opportunity to investigate more 
accurately the structure, properties, production mechanism as well as the inter- 
action of the f] meson with hadronic matter. The close-to-threshold total cross 
section measurements for the pp ppr] reaction [2-10], investigations on the 
differential cross sections for this reaction [8, 12, 39-41] and recently performed 
measurements of the analysing power for the pp — > pprj reaction [15, 16] made the 
theoretical analysis possible aiming in understanding the production process. 

3.1.1 Conclusions from the total cross section measure- 
ments 

Prom the total cross section measurements for the pp — > ppr) reaction it is inferred 
that amongst several possible scenarios the resonant current is the dominant one in the 
process of the rj meson production in proton-proton collisions. 

The majority of the theoretical models tries to elucidate the production of the 
rj meson within the framework of a one-boson-exchange formalism [11, 14,42-50], 
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where two interacting nucleons on the basis of a momentary energy violation 
exchange a virtual meson, which subsequently, under the interaction with nucleons 
turns into the 1] meson. Some of the possible mechanisms which may lead to the 
rj meson creation in nucleon-nucleon collisions are pointed out in Figure 13.11 and 
will be described in the further part of this section. What is characteristic to this 
kind of approach to the meson production is that the models do not enter into 
the quark-gluonic structure of the meson and nucleons, but rather introduce the 
phenomenological parameters like for example the coupling constants for different 
channels, the scattering lengths, the effective ranges, and cut-offs. 

On the other hand there are also trials to explain the r] meson production 
mechanism on the basis of instanton models for QCD vacuum [51,52]. In these 
Q CD-oriented approaches the calculations are on the very elementary level, where 
the effective degrees of freedom are not mesons and baryons like it was in the case 
of the one boson exchange models but rather the constituent quarks and gluons. 
These models however are temporarily in the early stage of development, and so 
we will deal here with the phenomenological meson-exchange models solely. 



N N N Nn NNj^lN 




a) b) c) d) 

Figure 3.1: Possible mechanisms of the rj meson production in nucleon-nucleon col- 
hsions: (a) resonant currents, (b) nucleon currents, (c) direct production, (d) mesonic 
currents. 

In Figure [371] there are depicted some of the possible scenarios for the 1] meson 
production in the nucleon-nucleon collisions. Figure I3.1[ a shows the resonant 
current, where the rj meson is produced in a two-step process: in the first stage 
the exchange of one of the pseudoscalar or vector mesons excites the nucleon to the 
Sii(1535) resonance and subsequently in the second step this resonance decays into 
a proton-?7 pair. Apart from the Sii(1535) resonance, also the contributions from 
the other possible resonances are considered. For example in the model reported 
in the reference [11] also Pii(1440) and Di3(1520) states, excited by the exchange 
of vr, 1], p, and u mesons have been included, however, the contribution from the 
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Sii(1535) excitation by the tt and t] mesons has been found to be dominant. The 
Sii(1535) resonance seems to play an important role as an intermediate state since 
it has a large width, covering the threshold energy for the pp pprj reaction, and 
couples strongly to the proton-?] system with the branching ratio corresponding 
to 30-55% [53]. The mechanism of the excitation of the colliding proton to the 
Sii(1535) resonance remains an open issue. The authors of [14,44,47] have found 
the p meson exchange to play the dominant role in this excitation. It has been 
stated that the p meson exchange is particularly important to explain the shape of 
the angular distribution of the pp —>■ pprj reaction. On the other hand, the authors 
of article [11] showed that the excitation function of the total cross section for the 
pp pprj reaction as well as the angular distribution of the emitted rj meson 
can be equally well described by the resonant mechanism, where the Sii(1535) 
is excited by the vr meson. This is presented in Figures I3.2[ a and I3.2[ b showing 
the comparison of the theoretical description of the excitation function for the 
pp pprj reaction: in Figure I3.2[ a the pseudoscalar meson exchange model has 
been used, while in Figure 13.21 b there are presented predictions of the vector 
meson dominance model. 
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Figure 3.2: Description of the close-to-threshold total cross section data for the 
pp ppr] reaction by means of the pseudoscalar meson exchange model (a) and vector 
meson exchange model (b). Data are from references [3-7]. The dotted curves are 
the resonance current contributions, the dashed lines represent the nucleonic current 
contribution, while the dashed-dotted curves are the mesonic current contribution. The 
solid lines show the total contribution. The figure is adapted from [11]. 



As can be seen from these figures, both models give equally well descriptions of 
the experimental data. The underestimation of the data by the models in the very 
close-to-threshold region results mainly from neglecting the r^-proton final state 
interaction (FSI). Proton-proton FSI has been accounted for in both calculations. 
Also at the higher excess energies (Q>60 MeV) models underestimate the existing 
data, which may be due to the higher partial waves which were not taken into 
account. The main conclusion here is that the resonant current (dotted curves). 
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i.e. the mechanism shown in Figure IXTl a. plays the dominant role in the t] meson 
creation. However, as far as only the excitation function of the total cross section 
is being considered, we are not able to state about the exchanged mesons in the 
process of the t] meson production in proton-proton collisions. 
Contribution of the other process that is depicted in Figure 13.21 is the nucleonic 
current (dashed curves) which we deal with when one of the nucleons emits a 
meson, which subsequently mixes its quark-gluonic structure with the structure 
of the other nucleon, and in the final state we again finish with a proton along 
with the r] meson (see diagram 13. li b). Naturally, the final state proton may have 
different properties - momentum and/or angular momentum, etc. - than the 
initial state proton. 

The last process taken into account in the abovementioned models is the mesonic 
current (dash-dotted curves, schematically presented in diagram 13. li d), when the 
rj meson is being produced via the exchange and fusion of two mesons emitted 
instantly from interacting nucleons. 

As has been mentioned previously, in this model the contribution of the two 
latter currents is rather negligible with comparison to the contribution of the 
resonant current (see Figure [3^ . 

Last possible scenario presented in Figure 13. li e - the direct rj production - 
although not taken into account in the calculations of [11] has been considered 
in the model of reference [50]. It has been found that the contribution from this 
process is of minor importance, being a one order of magnitude weaker than the 
mesonic current contribution. 



3.1.2 Differential cross sections 

Results of the measurements of the differential cross section for the pp pprj reaction 
are reported. 

Figure 13.31 shows one-dimensional differential spectra of the cross sections for 
the pp pprj process at Q = 36 MeV, namely the angular distributions of the 
emitted r] meson in the centre-of-mass (CM) frame. Left panel presents the theo- 
retical predictions of the pseudoscalar meson exchange model confronted with the 
data. The description of the curves shown in this figure is the same as that of 
Figure [3^ The data have the tendency to bend down in the range of extreme val- 
ues of the cos(6'). On the other hand, the contribution from the resonant current 
of the pseudoscalar meson exchange model (see Fig. I3.3[ a) has the opposite ten- 
dency, i.e. it bends upwards. However, the interference with the remaining two 
currents - nucleonic and mesonic - bends the overall contribution downwards. 
The behavior of the experimental data, as can be seen in Fig. I3.3[ a. is not quite 
well reproduced by the resultant interference. 
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Figure 3.3: Description of the angular distribution of the produced r] meson in the 
CM system in the pp ppr] reaction at Q = 36 MeV within the pseudoscalar meson 
exchange model (a) and vector meson exchange model (b). Data are from reference [39]. 
The meaning of the curves is the same as in Figure [3?2l The figure is adapted from [11]. 



In the case when the main contribution to the excitation of the Sn resonance 
originates from the p meson exchange (Figure 13.31 b) the experimental data are 
slightly better reproduced, especially at the edges of the cos(^^) range. This might 
indicate that the vector meson exchange mechanism was better in explaining the 
shape of the angular distributions for the pp — > ppr] process. 
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Figure 3.4: Angular distribution for the t] meson emission in the CM system for 
the pp — > ppr] reaction (a) as measured by the COSY/TOF collaboration [12] at 
Q = 41 MeV, (b) as obtained by the COSY-11 collaboration [8] at Q = 15.5 MeV 
(fuh circles) and by the COSY/TOF collaboration [12] at Q = 15 MeV (open circles). 
The solid line is the fit of the Legendre polynomials. Figures are adapted from [8]. 
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Here, a few words of explanation are required. At the time when the work 
of [11] has been released the only measurements of the differential cross section 
for the pp — >■ pprj reaction available were the early results of H. Calen et al. [39]. 
Unfortunately, the detector used in these measurements - the PROMICE/WASA 
detector [39] - had relatively poor angular acceptance in the side-range of 9 angle 
(i.e. the cos(^) around value). Later on, the measurements have been repeated 
by means of the COSY/TOF experimental setup [12] at the excess energy of 
Q = 41 MeV and show, within the error bars, the isotropic behavior over the 
whole range of cos^ - see Figure [3^ a. 

Therefore, as can be seen from comparison of Figures [3.31 and [3^ a. the theo- 
retical predictions suitable for description of the fiat differential cross sections are 
either those of the total contribution of the pseudoscalar meson exchange model 
or those of the resonance current contribution of the vector exchange model. 

The newer data on the angular distribution for the pp — > ppr] reaction at 
Q = 15.5 MeV [8] and Q = 15 MeV [12], both presented in Figure [331b, show the 
fiat distribution of the r] meson emission angle in the CM frame. The discussed 
angular distributions suggest that in the close-to-threshold region, i.e. up to the 
excess energy of about Q = 50 MeV, the rj meson is predominantly produced in 
the s wave and the higher partial waves are rather suppressed. 

3.1.3 Isospin dependence of the hadronic 77 meson produc- 
tion 

It is shown that the measurements of the isospin dependence of the total cross section 
for the A^A^ — > NNij reaction help to diminish the number of possible mesons in the 
excitation of the Sii(1535) resonance. 

The conclusion from the analysis of the results on the total and differential 
cross section for the pp — > pprj reaction presented in Sections 13.1.11 and 13 . 1 .21 is that 
the production of the 77 meson in proton-proton collisions proceeds predominantly 
via the excitation of one of the colliding protons to the resonant state Sii(1535). 
However, the relative contributions from the tt, rj, p, and u meson exchanges still 
remain to be determined. 

For the NN NNt] reaction there are two independent total cross sections: 
(Jo and (Ti, corresponding to the total isospin of two initial state nucleons 1 = 
and 1 = 1, respectively. The pp — > pprj reaction corresponds to the pure cxi cross 
section. The pn —>■ pnr] reaction is a mixture of I = and 1 = 1 initial states, as 
the proton's and neutron's isospins may sum up to the total isospin or 1. The 
theoretical considerations involving the Clebsch-Gordan coefficients yield: 

a{pn ^ pnr]) = ^{(ro + (Ti). (3.1) 

The total cross section for the quasi-free pn —>■ prirj reaction has been de- 
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termined in the excess energy range between Q = 16 and 109 MeV by the 
WASA/PROMICE collaboration [13]. The ratio R,, = a{pn pnri)/a{pp pprj) 
has been found to equal about 6.5 in the quoted excess energy range, as depicted 
in Figure [3751 
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Figure 3.5: The ratio R,, = a{pn pnrf)/a{pp pprf) of tiie total cross section 
for the quasi-free pn — > pnr] to the total cross section for the pp — > pprj reaction as the 
function of the excess energy. The original data are from reference [13]. 



From relation |3.H the fact that a{pp pprj) = cxi, and the definition of the 
ratio R^ = a{pn pnri)/a{pp pprj) we get: 

ao = {2R^ - l)(Ti. (3.2) 

Therefore, the measured ratio i?^ = 6.5 implies that the production of the ?7 meson 
with the total isospin 1 = exceeds the production with the isospin 1 = 1 by a 
factor of 12. 

Using the short range approximation, where the S wave cross sections are given 
by [13,54] 

(j{pp ppr\) = ai{pn pnr]) = C\U + tn + tp + Cn^/=i(0)P, (3.3) 

and 

aoipn ^ pnri) = C\ - 3U + tr, + 3tp - C| V/=o(0)|', (3.4) 

where stands for the strength of the different meson exchanges. The relative 
phases of the meson exchanges were taken to be real. The factor C, which includes 
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the phase space and particles interaction, was assumed to be the same for both 
isospins. The squared ratio between 1 = and 1 = 1 wave functions iplr) at the 
zero distance between the particles was calculated to be approximately constant 
and equals 0.8 [54]. Therefore the ratio 



which is a relatively large value, suggesting the dominance of the isovector meson 
exchange (vr and p) over the isoscalar meson exchange for the t] meson hadronic 
production. 



3.2 How do the analysing power measurements 
enable a distinction between the possible 
scenarios? 



The hitherto performed measurements of the analysing power for the pp pprj reaction 
are quoted. The available model predictions for the Ay at different excess energies as 
calculated based on different meson exchange models are shown. It is pointed out 
how the measurements of the analysing power enable to infer about the production 
mechanism of the rj meson in hadronic collisions. 

What has been described previously implicates that more limitations have 
to be added to the models in order to extract the information which meson - 
pseudoscalar tt meson or vector p meson - plays the most important role in the 
excitation of the intermediating Sii(1535) resonance. 

One solution would be verification of the polarisation observables - like the 
analysing power or spin correlation functions - given by different models. At 
present there exist two models that predict the energy dependence of the beam 
analysing power for the^ — >■ pprj reaction [11,14]. There are significant differences 
between the predictions based on these models visible in Figure 13. 6t e.g. the 
relative sign and the magnitude of the analysing power. Measurements of this 
observable might therefore serve in establishing the valid mechanism of the 7] 
meson production. 
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Figure 3.6: Predictions of the analysing power for the pp — > ppr] reaction as a function 
of the cosine of centre-of-mass polar angle of ry - cos{9j^) - at Q = 10 MeV (a) and 
Q = 36 MeV (b). Full lines are the predictions of the pseudoscalar meson exchange 
model [11] whereas the dotted lines represent the results of the calculations based on 
the vector meson exchange [14]. 

According to the model predictions of reference [14] the analysing power for 
the pp —>■ ppr] reaction in the close-to-threshold region should have the following 
form: 

A,(g,^,) = A7^'(Q)sin(2^,), (3.6) 

where 6'^ is the angle of the t] meson emission in the centre-of-mass frame. Equa- 
tion [32] has been obtained under the assumption of only p meson exchange in the 
excitation mechanism of the Sii(1535) resonance. The amplitude A^"'^ in refer- 
ence [14] is parameterized by means of the dimensionless parameter rj, namely: 

^rnax ^ _i.5[(_3 ± 0.5)7/^ + (16 ± 4)r/^], (3.7) 

where 

I 4M 

In Equation 13.81 m = 547.75 ± 0.12 MeV [53] denotes the rj meson mass, while 
M = 938.27203 ± 0.00008 MeV is the proton mass [53]. The shape in the form 
of sin(2^^) was obtained from the parameterization of the 'yp — >■ pr] differential 
cross sections [55,56] with the Legendre polynomials, and by the assumption that 
vector-meson-induced reaction amplitudes may be obtained from the amplitudes 
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for photoproduction, using the proper Jacobians [14]. The predicted values for 
jmax function of the excess energy Q are depicted in Figure 13. 7[ 



max,vec 




Q [ MeV ] 

Figure 3.7: Predictions for the value of A™*^^ for the pp — > pprj reaction as a function 
of excess energy. In the estimations, the vector meson exchange of reference [14] has 
been applied. 



The characteristic feature here is that A™"^ is peaked at Q = 10 MeV, where 
COSY-11 has the largest acceptance. This feature has been one of the main 
motivations to perform the measurement of the analysing power at this particular 
excess energy value. At Q of about 36 MeV A™"^ ^ 0, and we also decided to 
investigate this region of the excess energy. Note that at this value of excess energy 
the A^"-^ changes its sign and becomes negative for the higher excess energies in 
the close-to-threshold region. 

The model independent considerations of the spin observables can be found 
in reference [57]. The details of the pseudoscalar model assumptions and its 
prediction concerning the analysing power for the pp pprj reaction has been 
described in [11,58]. 

The first test measurement of the analysing power for the pp pprj reaction 
at the excess energy of Q = 40 MeV has been performed by the COSY-11 collab- 
oration in the year 2001. The method of the analysis and the results have been 
reported in [16,17] and are depicted in Figure 1X51 
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Figure 3.8: (a) Analysing power for the pp — > pprj reaction at Q = 40 MeV. Meaning 
of the curves is explained in the text, (b) Comparison of experimental data of the 
analysing power for the pp — > ppr] reaction at Q = 40 MeV with the predictions of 
pseudoscalar (solid line) and vector meson exchange models (dashed line). 

Figure [3^ a shows the comparison of the experimental results with the theo- 
retical predictions [59] of the individual partial waves contributions in the rj meson 
production. Production of the rj meson in the s wave solely would force analysing 
power to vanish [59], and therefore any non-zero result must involve the produc- 
tion of the 1] meson in higher partial wave. The dotted curve corresponds to the 
case, where the t] meson is created in the s + p waves, while the dashed curve 
represents the t] meson production in the s + p + d wave. The solid line is the full 
result of the model, where also the higher partial waves of the r] meson production 
are taken into account. Unfortunately the accuracy of the data does not allow to 
state firmly whether the rj meson is produced in the s or s + p wave and therefore 
to cope with this problem further investigations were required. One should also 
mention that the data of [16] need to be reanalyzed in order to get rid of the 
occasional error, however as has been checked in the analysis described in this 
dissertation, the corrected results should not differ from the results presented in 
Figure 13.81 of more than one standard deviation. 

Figure I3.8[ b shows the comparison of the data with theoretical predictions of 
the pseudoscalar meson exchange model (solid line) and vector meson exchange 
model (dashed curve). Both model predictions lie within around 2a distance from 
experimental data, and therefore at the level of accuracy obtained during the first 
measurement of the analysing power we were not able to distinguish between two 
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different hypotheses of the rj production. Further investigations were necessary 
and constitute the subject of this dissertation. 



Chapter 4 
Experiment 



The laboratory aparatus which enabled the experiments presented in this disser- 
tation comprises the storage ring COSY, which provides the accelerated polarised 
proton beams, and the detection setup COSY-11 used to register and identify the 
products of reaction. Both facilities will be described in the following sections. 

The COSY-11 detection setup has been described in details in many previous 
publications, therefore here we will present it only very briefly emphasizing the 
detectors relevant for measurements described in this thesis and aspects connected 
with the beam polarisation. 

4.1 Cooler Synchrotron COSY 
4.1.1 General properties 

The rough description of the cooler synchrotron COSY is given. 

The cooler synchrotron COSY [21,22] - a storage ring designed to accelerate 
the polarised and unpolarised beams of protons and deuterons - is operated by 
the Institute of Nuclear Physics (IKP) in the Research Center Jiilich in Germany. 
This device, sketched in Figure 14. 1^ consists of the ion source that provides the 
polarised or unpolarised ions of H~ and D~ to be further preaccelerated in the 
low energy cyclotron JULIC, which accelerates these beams up to the energies of 
45 MeV for H~ ions and 75 MeV for D~ [60]. After the preacceleration, the beam 
is extracted and guided through the 100 m long beamline to be injected into the 
storage ring COSY, that further accelerates the beams to the demanded momenta 
in the range from 0.3 GeV/c up to 3.7 GeV/c. 

The synchrotron consist of two arcs, containing in total 24 dipole magnets, 
and two straight sections, each of about 40 m length. The total length of the 
COSY ring is 183.4 m. 

The synchrotron COSY is equipped with two systems of beam cooling: the 
electron and stochastic cooling [70] . In order to increase the intensity of the beams 
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Figure 4.1: COoler SYnchrotron COSY - tiie floorplan. Several internal (ANKE [61], 
COSY-11 [62], PISA [63]) as well as external experiments (GEM [64], HIRES [65], 
MOMO [66], TOF [67,68]) are making use of the COSY beams. Recently the Att 
detection setup WASA [69] has been installed. 
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of polarised protons the electron cooler is being used. Using the electron cooler and 
the stacking technique it was possible to store 10^° polarised protons during the 
runs concerned in this thesis. The stochastic cooling for COSY [71,72] is designed 
to reduce the emmitance of the proton beams in the momentum range between 
1.5 and 3.7 GeV/c. The system consists of one pickup tank of 4 m length and a 
kicker of 2 m length for each the vertical and horizontal planes. The stochastic 
cooling is generally used for internal target experiments in order to achieve the 
equilibrium conditions between target heating and stochastic cooling. For more 
details on the working principle of the COSY electron and stochastic cooler the 
reader is referred to [70,72]. 

4.1.2 Production of the polarised proton beam 

The source of the polarised proton beam is described along with the method of polarised 
beam production at the cooler synchrotron COSY. 

The source of the polarised proton beam [73] is of a Colliding-Beam Source 
type, that provides the polarised H~ ions in a direct charge-exchange process of 
colliding beams: 

+ ^H- + Cs+. (4.1) 

Here, the neutral nuclear polarised beam of hydrogen produced in the atomic- 
beam source meets a fast neutral Cs" beam. The cross section for the process 14.11 
is large because the binding energy of the electron in Cs (3.9 eV) is close to the 
binding energy of the electron in hydrogen's 2S state (3.4 eV). Additionally the 

atoms are rather slow, which increases the probability of the charge-exchange 
reaction. For injection the nuclear polarised anions are used rather than H"*", 
as it has been experimentally proven that stripping off the electrons from H~ 
during injection into the COSY is about an order of magnitude more efficient 
than stacking injection of protons using a bunched beam [73]. 

The polarised ion source, presented in Figure IT2l consists of dissociator, two 
groups of sextupole magnets, the radio frequency transitions, the solenoid magnet, 
cesium ionizer, colliding zone, the 90° deflector, and a Wienfilter. In the dissocia- 
tor the H2 molecules are dissociated into atoms, which subsequently through the 
cooled nozzle reach the sextupole magnets region. The first sextupole magnet pro- 
duces the electron state polarisation of the atoms, and subsequently atoms with 
the electron spin state mj = —1/2 are defocused and only atoms with mj=+l/2 
survive. The second sextupole magnet acts as the lense to focus the hydrogen 
atoms into the radio frequency transition units, which produce the nuclear polar- 
isation. Afterwards the polarised atoms collide with a cesium beam moving 
in the opposite direction. In the charge-exchange region where the reactions 14.11 
occur the nuclear polarisation of hydrogen is preserved by the strong longitudinal 
magnetic field at the order of 1.7 kC. For this value of magnetic field the polarisa- 
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Figure 4.2: The polarised ion source for COSY. Figure is adapted from [73]. 

tion of protons at the level of circa 95% can be expected [73]. The polarised H~ 
anions are subsequently deflected by 90° in the magnetic deflector and guided into 
the extraction system, passing on the way the Wienfilter, which may be rotated 
around the beam axis, and therefore can select the required spin orientation of the 
H~ ions, and also separate the anions from the electrons and other background 
particles. 

Subsequently the H~ beam is injected into the cyclotron JULIC at the energy 
of 4.5 keV. The cyclotron preaccelerates them up to energies of 45 MeV. After 
that, the beam of H~ anions is guided through the beam lines and reaches the 
stripping injector [74], where the H~ anions are stripped off two electrons and 
furthermore injected into the COSY-Ring as polarised protons with the selected 
polarisation orientation. 

With this method the intensity of circa 10^° stored polarised protons with a 
degree of polarisation over 65% have been obtained [75] during the experiments 
reported in this work0. 

4.1.3 Acceleration of the polarised proton beam 

Acceleration of tiie polarised proton beams is described. Superimposed are the meth- 
ods of overcoming the intrinsic and imperfection resonances. The table of all possible 

^Nowadays it is possible to achieve over 90% polarisation for the internal experiments and 
circa 80% for the external ones [76]. 
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depolarising resonances at the COSY ring is presented. 

The acceleration of the polarised beams is challenging due to the number of de- 
polarising resonances that occur during this process, depending on the demanded 
momentum of the beam. 

The polarisation vector of the beam of protons stored in the ring precesses 
around the magnetic fields encountered along the particles' orbill^. In the ideal 
ring, where the dipole magnetic fields are vertical, the number of spin precessions 
during one turn (so called spin tune) is given by [77]: 



where 7 is the Lorentz factor, and G denotes the gyromagnetic anomaly - for 
protons G = 1.79285. Depolarising spin resonances arise when the frequency 
of the precession of the polarisation vector is such that upon each revolution 
the difference between the phases of the polarisation vector and the depolarising 
magnetic field vector is the same at the point where these two meet. At the COSY 
ring there are two main types of depolarising resonances [75]: 

• imperfection resonances, caused by magnetic field errors, arising if the 
magnets are slightly misaligned or if there are vertical orbit distortions. 
These resonances occur at the beam momenta for which 76* has an integer 
value; 

• intrinsic resonances, caused by the radial fields due to vertical focusing of 
the beam, occurring when •yG ± {Qy — 2) equals an integer, where Qy is the 
betatron vertical tune, i.e. the number of vertical oscillations of the beam 
per one revolution. 

The resonances connected with the horizontal betatron tune Qx, occurring when 
'jG ± {Qx — 2) equals an integer are not important at COSMj for the vertically 
polarised beams, because the horizontal betatron motion is driven by the vertical 
magnetic fields of the quadrupoles. Also we will not consider here the induced 
resonances, arising due to the longitudinal oscillating magnetic fields. 

Table 14.11 [78] shows the chart of all depolarising resonances and the corre- 
sponding proton beam momenta Pbeam in the momentum range of COSY. Reso- 
nance strengths and ratios of preserved polarisation Pf/Pi are shown. In the 
simulations [78] an energy gain per turn equal to 0.7 keV has been assumed. Num- 
ber of intrinsic resonances depends on the superperiodicity S, the parameter which 
describes the machine setting of the quadrupole magnets. If all magnets are oper- 
ated with the same quadrupole settings, the superperiodicity S = 6 and only one 
intrinsic resonance occurs. However, this setting does not allow to accelerate the 

^In the case when the precession repeats on each revolution the particles are on the so called 
spin closed orbit. 

■^As long as the coupling between the horizontal and vertical betatron motion is not present. 
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beam up to the maximum beam momentum. For this purpose a special setting of 
the magnets is required [78] with S = 2, thus four additional intrinsic resonances 
appear, all listed in Table 14.11 The occurance of the imperfection resonances does 
not depend upon the superperiodicity. 
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-1.00 
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2781.2 


0.53 


-0.74 
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2996.4 


0.46 


-0.58 


INT 


10-Q^ 


2 


3208.9 


0.25 


0.25 



Table 4.1: Imperfection (IMP) and intrinsic resonances (INT) at the synchrotron 
COSY. 

The ratio of the preserved polarisation Pj/Pi depends on the strength of a 
resonance (see Table I^TTi) . Each time the beam is crossing the imperfection reso- 
nances, the polarisation may be completely lost if the strength of the resonance 
is not sufficient to flip the spin of all beam particles. At some e^- value, for a 
certain imperfection resonance and certain momentum spread of the beam Ap/p, 
the ratio Pf/Pi is equal to —1, which means that crossing the resonance with a 
proper strength one deals with the spin flip, with no polarisation loss. Simula- 
tions reported in reference [78] showed, and the experiments confirmed, that the 
excitation of the vertical orbit by 1 mrad, using the horizontal correcting dipoles, 
is sufficient to adiabatically flip the spin at all imperfection resonances. 

Another method to overcome the imperfection resonances involves the use of 
spin rotators - so called Siberian snakes. At the synchrotron COSY the solenoids 
of the electron cooler, acting as the partial Siberian snake, are able to rotate 
the spin around the longitudinal axis by circa 8° at the maximum momentum 
of COSY [78]. It has been experimentally proven that a rotation angle of less 
than 1° leads to a total spin flip with no polarisation loss at all flve imperfection 
resonances. 

In order to avoid the intrinsic resonances the technique of tune jumping is 
applied. A sudden tune jump increases signiflcantly the crossing speed of the 
resonance and therefore the beam "jumps over" the intrinsic resonance with less 
than 5% polarisation loss at the strongest intrinsic resonances and no more than 
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1% loss at all other intrinsic resonances [78]. 

The experiments discussed in this dissertation have been performed at the 
beam momenta of pteam = 2.010 and 2.085 GeV/c. Thus, according to Table WA\ 
three imperfection and two intrinsic resonances had to be crossed. Finally, as we 
shall see in Section 15.41 the obtained polarisation degree was equal to about 68% 
and 66% at the lower and higher beam momentum, respectively. 

4.2 COSY-11 facility 

The detectors needed for identification of the pp ppr] and pp pp reactions are 
described. The trigger conditions for both reactions are given. 

The experiments described in this dissertation have been performed by means 
of the COSY-11 detector setup [18] presented in Figure 14.31 COSY-11 is an 
internal experiment at the cooler synchrotron and storage ring COSY, designed 
to study the production processes, the structure, and interactions of the mesons in 
the 1 GeV mass range. The facility makes use of a regular COSY dipole magnet, 
acting as a magnetic momentum spectrometer for charged reaction products. 

A vertically polarised proton beam has been scattered on the H2 molecules 
from an internal cluster target [79, 80] installed in front of the COSY magnet. 
Reaction products possess lower momenta than the beam protons, therefore these 
are bent more in the magnetic field of the dipole. Positively charged ejectiles leave 
the scattering chamber [18] through the thin exit windoAMj reaching the detection 
system operating under the atmospheric pressure. 

Trajectories of the positively charged protons, which are bent to the inside 
part of the ring, are measured by means of two planar drift chambers [18,81] Dl 
and D2, presented schematically in Figure 14. 3[ These drift chambers are spaced 
by 70 cm and contain together 6 + 8 = 14 detection planes. The active area of 
the chambers is 1680 x 433 mm^. The Dl chamber contains 6 detection planes; 
two planes with vertical wires, two with wires inclined by +31°, and last two with 
wires inclined by —31°. The D2 drift chamber has additional two planes in the 
back part with vertical wire orientation. The wires of the consecutive planes of 
each pair are shifted by half of the cell width (20 mm) in order to resolve the 
left-right ambiguity of passing through particles. Drift chambers operate with 
a gas mixture of 50% argon and 50% ethane at slightly more than atmospheric 
pressure. The particle trajectories are reconstructed using the computer code 
MEDUZA [82]. The code allows for reconstruction of multiple track events. 

^This exit window ( 1870 x 76 mm^) is made of a 30fim layer of aluminum foil, and two crossed 
unidirectional sheets of carbon fibers soaked in epoxy resin of f 50/im width, acting as an outward 
carrier material. The choice of materials with low nuclear charge reduces straggling in the exit 
window to values below the resolution of the detection system [f8]. 
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FigTirc 4.3: Schematic top view of the COSY-11 detection setup [18]. Only detectors 
needed for measurements of the pp — > pprj and pp ^ pp reactions are shown. SI, S3, 
and S4 denote the scintillator detectors; Dl and D2 stand for the drift chambers; Si,„on 
are the silicon detectors for measurements of the elastic scattering. Figure is adapted 
from [23]. 

The particle trajectories are traced through the magnetic field of the dipole 
back to the interaction point. In this way the momenta of the particles can be 
determined with precision of 4 MeV/c (standard deviation) [8]. 

Leaving the drift chambers, the positively charged reaction products pass the 
scintillator hodoscope SI [18], consisting of 16 scintillation modules with dimen- 
sions of 450 X 100 X 4 mm^ made of Bicron BC 404 [83]. The scintillators are 
arranged vertically and read out at both sides by photomultipliers. The time 
resolution of this hodoscope in the range between 160 and 220 ps has been mea- 
sured [84], depending on the hit position. This detector dehvers the start signal 
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for the time-of-flight measurement. 

The stop signal of the time-of-fiight measurement is generated in the S3 scintil- 
lator wall [18] (made of a 220 x 100 x 5 cm^ scintillator Bicron BC 404), situated 
at a distance of circa 9 m from SI. The scintillator wall generates a light signal 
upon a charged particle crossing its volume. This signal is read out by a matrix of 
217 photomultipliers. The scintillator wall and the photomultipliers are separated 
with 4 cm air gap. The centre of gravity of the pulse height from individual pho- 
tomultipliers is calculated in order to resolve the hit position of a particle. The 
restriction of few photomultipliers responding allows to separate two or more hit 
positions in the scintillator wall. 

Time-of-flight measurement together with a reconstructed momentum of a 
particle yields a possibility of indentification of the particle by calculating its 
mass. Hence, the four momenta of the outgoing particles can be determined. 
This information, in connection with the known four momenta of initial state 
particles allows to calculate the missing mass of the undetected particle or system 
of particles. The missing mass method will be described in Section 15.1.51 

The trigger condition for registering the pp — > pprj reaction is: 

) A S3^>2, (4.3) 

where fi denotes the multiplicity of segments in the SI detector, and also the mul- 
tiplicity of the photomultipliers that have fired in the S3 detector. Subscript high 
stands for the high energy deposition in the SI scintillator module, corresponding 
to the case when two or more particles cross a single module. 

Similarly, for registering the pp — > pp reaction a trigger condition of the type 

Tp-p^pp = 51^=1 A SA (4.4) 

has been set, where S4 is a scintillator detector used for registering the recoil pro- 
tons from the elastic scattering. Protons passing through this detector reach the 
granulated silicon-pad detector Simon consisting of 144 silicon pads with dimen- 
sions of 22x4.5x0.28 mm^ arranged in three layers one above the other. Two of 
the layers are located in the front part of the detector and one in the back. 

For completeness we would like to mention that there are other detectors build- 
ing up the COSY-11 facility like the silicon pad detector [85,86] for the spectator 
protons from the quasi-free pd — > ppnX reactions [87-89], the neutron detec- 
tor [85,90-92], the Cherenkov detector [93], the C-shaped hexagonal chamber [94] 
and the auxiliary scintillator detectors [18,23]. However, as these detectors were 
not used during the experiments reported in this dissertation we omitted their 
description herein. The interested reader is referred to the publications quoted 
above. 
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5.1 Calibration of the detection system 
5.1.1 Time-space calibration of the drift chambers 

The method of time-space function derivation for the drift chambers is described. 

The drift chambers described in Section 14.21 operate with a gas mixture of 
50% argon and 50% ethane at shghtly more than atmospheric pressure. Upon 
crossing of a charged particle through the gas-filled volume the electron clusters 
are generated and move towards the anode wires. The drift time of the electron 
clusters to the anode wire is a measure of the minimum distance between the 
sense wire and the trajectory of the particle. This relation is called a time-space 
calibration of the drift chamber and has to be determined from the experimental 
data. 

The drift velocity varies with the atmospheric conditions [95]. Therefore, in 
order to perform the time-space calibration the experimental data have been di- 
vided into 32 groups for each energy, each group corresponding to about 6 hours 
of measurement. The time-space calibration [23] has been performed iteratively 
starting with an approximate time-space function d{t). From this function and 
from measured drift times the points where the particle crossed the middle of the 
cells have been found for each detection plane and the straight line corresponding 
to the particle's trajectory has been fitted. For such obtained trajectory of the 
particle for an i-th event the distance d[ to the sense wire has been calculated and 
the average of the differences 



obtained from the mentioned sample of the experimental data have been used 
in order to correct the time-space function. In the Equation 15.11 dt denotes the 
distance between the particle trajectory and the sense wire for an i-th event. 
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n 



38 



Data Analysis 



reconstructed using tlie d(t) calibration function. Subsequently, new time-space 
functions d{t) + Ad{t) have been used for further iteration. The procedure has to 
be repeated until the Ad corrections becomes negligible. 




drift time + offset [ ns ] 

Figure 5.1: Distance from the sense wire as a function of the drift time for an arbitrarily 
chosen ceh of DCl - so cahed time-space cahbration of the drift chambers. 

Figure [STTl a presents the time-space calibration of an arbitrarily chosen sense 
wire of the DCl. The linearity of the distance from the sense wire as a function of 
the drift time may be seen in the range between circa 100 and 500 ns. The value 
of 100 ns corresponds to the time offset introduced by the electronics. 

5.1.2 Time-of-flight calibration 

Time calibration of the scintillator detectors used for time-of-flight measurement is 
presented. 

As it was mentioned in Section 14.21 the SI and S3 counters serve as a start 
and stop detectors, respectively, for the time-of-fiight measurement. As the SI 
scintillator consists of 16 modules, in which signals are read from both sides by 
photomultipliers and the signals from S3 scintillator are read by a matrix of 217 
photomultipliers, in order to obtain the credible time-of-fiight information, the 
relative time offsets of the electronics for each photomultiplier have to be deter- 
mined. 
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For the SI time cahbration, we used the events where particles were crossing 
adjacent modules. In this case we can assume that in both modules the signal is 
generated at the same time. The TDC value of the single photomultiplier of the 
SI detector can be expressed as follows [23]: 

TDCsl = t.l +ty+ t^f\PM) + tt^''\PM) - ttrigger, (5.2) 

with tsi denoting the real time of the signal generation in the SI detector and 
ty standing for the time that the signal needs for passing the distance between 
the hit position in the scintillator module and the edge of the scintillator. The 
t^°''^(PM) is the time walk effect correctioE0, t°{'^'*^*(PM) denotes the time offset 
of the electronics for a given photomultiplier, and ttrigger is the time when the 
trigger pulse started the readout of TDC modules. 




Figure 5.2: (a) Determined TDC time difference between arbitrarily chosen 6^^ and 
gth mQfiyig of the SI scintillator detector after the calibration, (b) Difference between 
time-of-flight measured on the S1-S3 path from the signals registered in SI and S3 
detectors, and the time-of-flight between SI and S3 calculated from the reconstructed 
momentum of the particle. 

After the drift chamber calibration, we have chosen the events with only one 
track reconstructed and a signal in two neighboring modules of the SI detector. 
The time offset for the first module in SI (the one standing most to the inner part 
of the ring) has been set arbitrarily. Then, the time difference: t^i 2*"^* — TDC5i,2 
- TDCs'i,! has been taken as a correction for the time offset between the second 

^By the time walk effect we understand the variation of the registered TDC time with the 
amplitude of the signal typical for leading edge discrimination [96] . This effect can be corrected 
by adding an offset t™°"^(PM)= ^==, with j4_DC denoting the signal charge value [23]. 
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and the first module of the SI detector. The same procedure has been repeated 
for all modules. Figure [5^ a depicts the TDC time difference between arbitrarily 
chosen 6**^ and 5^^ module of the SI counter upon single particle crossing through 
the overlap region, after the correction for the time offset of the electronics and 
photomultipliers. One can see that this time difference is peaked around the zero 
value, which confirms the correctness of the time calibration of the SI scintillator. 

Having adjusted the time offsets of the individual photomultipliers in the SI 
detector and hence having established the common start for the time-of-fiight 
measurement we need to adjust the time offsets for each of the 217 photomultipli- 
ers in the S3 detector in order to set up the common stop. Similarly as it was in 
the case of the SI detector, the TDC value for the individual electronic channel, 
corresponding to the single photomultiplier in the S3 detector reads: 

TDCsS = ts3 + tpos + rj'\PM) + t%^''\PM) - ttrigger. (5.3) 

Here, the tpos denotes the time that light signals need to pass from the scintillation 
origin down to the photomultiplier photocathodes. 

In Equations 15.21 and 15.31 the ttrigger values are the same, and so the time-of- 
fiight value 

T0F{PM) = t,3{PM)-t,^ (5.4) 

does not depend on the triggering time. Hence, the only unknown quantities in the 
TOF calculation are the time offsets for the individual photomultipliers in the S3 
detector - t°l^^^^{PM). These values determined for each photomultiplier can be 
extracted from comparing the TOF values with the time-of-fiight calculated from 
the reconstructed momenta of the particles. The differences after the calibration 
for the photomultipliers no. 140-160 are presented in Figure I5.2[ b. In order to 
extract the t°l^'^'^^{PM) values the one-track events with an identified proton has 
been chosen. 

Having established all the time offsets of the photomultipliers in the SI and 
S3 scintillators, the time-of-fiight were calculated as TOF = tgs — tsi. One should 
mention that the t^i have been taken as the average of times from the upper and 
lower photomultipliers of the hit module, and t^s were the weighted mean times 
for the proper cluster of photomultipliers in the S3 detector that registered the 
light pulses. For more details on the SI and S3 calibration the reader is referred 
to [23]. 

5.1.3 Position of the drift chambers 

A method for the search of the optimal position of the drift chambers is presented. 

The time-space calibration of two drift chambers Dl and D2 has been de- 
scribed in Section I5.1.1[ It is very important to know the exact position of these 
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drift chambers in order to achieve the required resolution in the momentum recon- 
struction of outgoing protons. This position can be parameterized for example 
introducing the parallel shift of drift chambers, which in the following we will 
denote by Ax, and their inclination Aa. Both parameters are illustrated in Fig- 
ure 15.31 



DCl 



DC2 




a) b) 

Figure 5.3: Schematic illustration of the parameters Ax (a) and Aa (b) mentioned 
in the text. 

The parameters Ax and Aa influence the resolution of reconstructed four mo- 
menta of outgoing particles, and consequently influence also indirectly the missing 
mass resolution a {mm?). The studies of the missing mass resolution as a func- 
tion of Ax and Aa have been performed on the data set for the excess energy 
of Q = 10 MeV in order to optimize the position of drift chambers. We have 
varied these parameters in the range of [—1cm; 2cm] for Ax, and [—0.015°; 0.055°] 
for Aa, investigating the width of the missing mass spectra as a function of 
these parameters. The minimum value of the width of the t] peak in the missing 
mass squared histogram (which we will present in Section 15.1.51) has been found 
for Ax = 0.35 cm and Aa = 0.025° and equals a {mm^) = 0.00164 GeVVc^. 
This corresponds to the value of the width of rj missing mass peak which is 
a (mm) = 1.5 MeV/c^. 

For the analysis at the excess energy of Q = 36 MeV the same optimal param- 
eters have been chosen. 

It is important to note that the determined values of Ax and Aa are in a very 
good agreement with the result obtained in the preceeding COSY-11 analyses [17, 
25]. 

5.1.4 Beam-target relative settings 

It is shown that the relative position of the beam and target influence the reconstructed 
momenta of the particles. A search for the optimal setting is presented. 

Another parameter that influences the resolution of the reconstruction of four 
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momenta of protons are the dimensions of beam and target and also the relative 
position between them. In the procedure of momentum reconstruction it is as- 
sumed that the interaction vertex is a point where the infinitesimally thin target 
intersects with the absolutely thin beam of protons. However, in reality both beam 
and target are of finite dimensions. Therefore, the reactions may be initiated in 
some finite volume, size of which depends on the abovementioned parameters, 
introducing a spreading to the reconstructed four momenta of protons. In the 
analysis, following the studies reported in [27,97], we used the approach that 
the target may be described by a cylinder of radius r = 4.5 mm with uniformly 
distributed protons [79,80], and the density of the protons in the beam is given 
by the two dimensional gaussian distribution, as schematically illustrated in Fig- 
ure 15. 4[ The ax and cxy denote spreadings of the beam distribution in horizontal 
and vertical directions, respectively. Basing on the previous measurements [97], 
in our Monte Carlo studies we assumed these parameters to equal ax = 0.2 cm 
and ay = 0.5 cm. Distance between the centre of the target and the centre of 
proton beam distribution is denoted in Figure WM a, by hx- 



TARGET 




Figure 5.4: Schematic illustration of the beam and target dimensions and their relative 
position: (a) view from above, (b) side view. Figure is adapted from [27]. 



Let us denote by Ax the deviation between the center of the reaction region 
and the nominal position of the target. In order to determine the position of the 
center of the reaction region we used the elastically scattered protons, registered 
by means of the detectors shown in Figure [5751 Forward scattered protons from the 
elastic process are detected by the scintillator detector SI and the drift chambers 
Dl and D2. The recoil protons passing through the scintillator S4 reach the 
granulated position-sensitive silicon detector Simon- 
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Figure 5.5: Schematic view of the COSY-11 detection setup. Only detectors used for 
the measurement of the proton-proton elastic scattering are shown. Figure is adapted 
from [23]. 

As the proton scattered forward under the 62 angle - passes through the stack 
of drift chambers Dl and D2, its trajectory may be reconstructed and subsequently 
traced back to the reaction vertex, through the known map of the magnetic field 
inside the COSY dipole. This procedure yields the momenta of the particles pass- 
ing through the drift chambers. Further, the reconstructed momentum vector may 
be decomposed at the reaction point into two components: the transversal (p_l) 
and the longitudinal one (py ) , with respect to the momentum of the beam protons. 
The two body kinematics puts some constraints on these two components, namely 
the set of points should form an ellipse [98] in the momentum plane. The 

size of this ellipse depends on the y/s - the total energy in the centre-of-mass 
system. A part of one arm of this theoretical ellipse for the ^/s = 2.43 GeV, cor- 
responding to the proton beam momentum of Pbeam = 2.010 GeV/c, is depicted 
in Figures [5. 61 a and 15. 61 b as the solid line. 

The value of the reconstructed particle's momentum depends on the correct- 
ness of the assumption of the reaction vertex position. Shifting the reaction vertex 
upwards the beam line presented in Figure [5751 would yield in higher reconstructed 
momenta of the particles, as we decrease the curvature of the trajectory inside 
the dipole magnet. On the other hand, a parameter Ax with the negative values 
(shift of the reaction vertex downwards in Figure 15.51) would induce a higher cur- 
vature of the tracks inside the dipole magnet, hence lower momenta of the reaction 
products. The values of the components of the reconstructed momentum vector 
correspond to a point in the momentum plane p^ - Py. This point is lying inside 
or outside the expected theoretical ellipse, depending whether the reconstructed 
particle's momentum is lower or greater than the actual one. 
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Figure 5.6: Perpendicular versus parallel momentum component of the elastically scat- 
tered protons measured at the beam momentum of Pbeam — 

2.010 GeV/c for Ax = 

(a) and Ax = 0.7 cm (b). The solid line represents the theoretical kinematical ellipse 
for this particular value of the beam momentum. Enhancement of the event distribu- 
tion shows the elastically scattered protons over the constant background originating 
from the production reactions. Figures (c) and (d) show the projections of the event 
distribution along the theoretical ellipse for data from figures (a) and (b) respectively, 
(e) Distance between the theoretical ellipse and the centre of the points distribution 
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Following the studies of [97] we have been searching for the optimal position 
of the reaction vertex, by changing the Ax parameter in the range from to 
4 cm. In Figure 15.61 a a distribution of the reconstructed momenta of the elasti- 
cally scattered protons is presented, for the case where Ax was set to Ax = 0. 
The enhancement of the event distribution is assigned to the elastically scattered 
protons. The background is mainly due to many body reactions. As one can see, 
the reconstructed momenta of the elastically scattered protons are lying inside 
the theoretical kinematical ellipse, which means that the reconstruction yields 
too small values of the particle's momentum. This is also seen in Figure 15.61 c. 
where we presented the projection of the event distribution along the kinematical 
ellipse. One can notice the shift of the centre of the distribution outside the zero 
value. 

For the same data sample and with Ax = 0.7 cm the components of the 
reconstructed particle's momenta are lying on the kinematical ellipse, as can be 
seen in Figure I5.6[ b. Figure I5.6[ d shows the projection of the data set onto the 
kinematical ellipse for this case. Indeed, the event distribution is centered around 
the zero value. This means that for the data sample we considered here, the 
optimal value of the Ax parameter equals Ax = 0.7 cm. 

The same procedure has been differentially performed for the whole sample 
of data. The experimental data have been divided into 32 groups, corresponding 
to about 6 hours of measurement, and for each out of these 32 groups the value 
of the Ax parameter has been determined using a test, where was calcu- 
lated between the experimental event distribution and the theoretical kinematical 
ellipse. 

In Figure I5.6[ e we have presented the relation between the parameter Ax 
and the shift of the centre of the peak, originating from the projection of the 
experimental points along the kinematical ellipse. This dependence have been 
studied in the Ax range from to 4 cm, and it was found that it may be described 
fairly well by a polynomial function of second order, the one presented in the 
Figure I5.6[ e. One can see in this figure that the zeroth shift of the peak centre 
corresponds to the value of Ax ~ 0.7 cm. 

5.1.5 Missing mass technique 

The idea of the missing mass technique is sketched. The software presort of the events 
with two protons in the final state is given. We also present the missing mass spectra 
for both excess energies. 

As the ri meson is a non-charged short living particle, its direct registration in 
any of the nowadays known detectors is impossible due to the very low distance 
that is passed by this particle before it decay^. To our best knowledge, nowadays 

^The Particle Data Group [99] gives the value of the total width of the r\ meson F,, = 
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three methods are apphed in the worldwide experiments in order to identify the 
7] meson. First of them is the missing mass technique, which will be described in 
the following part of this section. Second method is the reconstruction of the four 
momentum of the rj meson by means of the identification of its decay products, e.g. 
27 rays. Third method is the combination of the two latter techniques, namely, the 
measurement of the 27 decay in coincidence with the directions (and/or energies) 
measurements of two final particles. This method, for example, will be used at 
the recently installed WASA-at-COSY experiment [69]. 

The idea of the missing mass technique involves the theorem of the four mo- 
mentum conservation. Let us consider the pp — > ppX reaction with the acceler- 
ated proton beam and the proton target with zero momentum value. Denoting 
by Pb = {Eb,Pb), Pt = (-E'tjO), Pi = (i?i,pi), and P2 = (-^2,^2) the four momenta 
of the beam, target, and two outgoing protons, respectively, one can calculate the 
missing mass (mx) of an undetected particle or system of particles in the exit 
channel according to the formula: 



After the adjustment of the drift chambers' parameters Ax and Aa and the 
parameter A^, and having performed the time calibration of the SI and S3 scintil- 
lators, as well as the calibration of the drift chambers (described in Sections 15.1. 11 - 
I5.1.4P one can reconstruct the momenta of the outgoing particles and also their 
velocities. Having the momenta (p) and velocities {(3) of the particles one can 
calculate their squared invariant masses applying the relativistic formula: 



Figure 15.71 shows the distribution of the invariant masses of two particles mea- 
sured in coincidence at the excess energy of Q = 10 MeV. Over the background 
arise four "isles" of the events with 7r+-7r+, 7r+-proton, proton- vr"*", and proton- 
proton identified in the final state. For the further analysis only events with two 
protons were taken into account. The software cut on the squares of the missing 
masses of particles is denoted with the solid lines. The dotted lines show the 
values of the tt"*" meson and proton squared masses. 



1.30 ± 0.07 keV, which yields the mean life-time of the ry meson equal to r,, = (5.10 ± 0.29) • 
10~^^ s, corresponding to cr,, = 1.53 • 10^^" m. 



ml = El -pI = (A + Pt-Pi- P2f = 
{Eb + Et- El- E2f - {pb + Pt-Pi- ■ 



(5.5) 



(5.6) 
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Figure 5.7: Squared masses of two simultaneously measured particles in the exit 
channel. Events from the measurement of the pp — > pprj reaction at the excess energy 
of Q = 10 MeV are presented. 



Knowing all the four momenta - P^, Pt, Pi, and P2 - and applying Equation l5.5l 
one can calculate the missing mass of X for each event separately. The spin- 
averaged missing mass spectra for the pp pprq reaction as measured using the 
COSY-11 detection setup are presented in Figure [531 Over the wide background, 
originating mainly from the multipionic production, clear rj peaks are visible for 
both measurements. There are around 3000 and 1500 rj events at Q = 10 and 
36 MeV, respectively, integrated over both spin orientations. 



1500 







PP^PPTl 

(Q = 10MeV) 




0.05 0.1 0.15 0.2 0.25 0.3 



missing mass^ [GeV^/c"] 



a) 



800 



o 

> 600 



10 

c 400 
> 



200 







PP^PPTl j' 




(Q = 36MeV) J 





0.26 



0.28 



0.3 



0.32 



missing mass^ [GeV^/c"] 



0.34 



b) 



Figure 5.8: (a) Spin-averaged spectrum of the square value of the missing mass for 
the pp — > ppr] reaction at the excess energy of Q = 10 MeV, as measured by means of 
the COSY-11 detector setup, (b) The same, but for the excess energy of Q = 36 MeV. 
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The differences in the shape of the background are due to the different trigger 
setting between two measurements, and also due to the fact that in both analyses 
different software cuts have been applied. It is also visible that the signal from 
the 7] meson production at Q = 10 MeV is closer to the kinematical limit than the 
signal at Q = 36 MeV. This is due to the fact that the kinematical limit m^ax 

rrimax = mn + Q, (5.7) 

equal^ m^ax= 0-311 MeVVc^ at Q = 10 MeV and it equals 
m^„^= 0.340 MeVVc^ at Q = 36 MeV. 



5.2 Identification of the spin up and down 
modes 

It is shown how we resolve the problem of the identification of the spin state of the 
polarised proton beam using the elastically scattered events. 

A few words of explanation are required concerning the method of spin up and 
down identification. In principle the polarisation state could have been determined 
from the known spin orientation in the polarised ion source and the number of 
imperfection resonances, mentioned in Section 14.1.31 to be crossed in order to 
accelerate the beam up to the required momentum. However, even without a 
precise knowledge of the beam optics during the acceleration process it is possible 
to determine the spin state of accelerated protons. 

To enable an offline assigment of the polarisation mode two scaler channels 
have been used, working in the sequence mode. Whenever there was a beam in 
the spin state A one of the scaler channels has been working while the other one 
remained idle, and vice- versa. 

In order to identify the spin state A with spin up or down mode the detector 
system of Figure 15.51 for registration of the elastic scattering in the accelerator 
plane has been used. 

For the elastically scattered protons at the beam momentum of pbeam = 
2.010 MeV/c, which reach the drift chambers Dl and D2 (see Figure [53!) the 
scattering angle 62 may vary from 35° to circa 80° in the centre-of-mass system. 
Here we deal only with scattering to the right with respect to the polarisation 
plane. For spin up mode, where the polarisation vector of the proton beam is 
pointing along the Oy"''^'^ axis defined in Section [2!2l the formula 12.131 - integrated 

•^In the calculations the meson mass = 547.75 MeV/c^ [53] has been used. 
^Which is restricted by the acceptance of COSY-11 setup for the pp pp reaction - see also 
Figure EH 
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over TUpp, TUpj^, and ifj angle - reads: 

ai92,P) = aoi92){l-PAy{92)) for cos ^2 G (35°, 80°) , ^ 180°; (5.8) 
while for spin down it yields: 

ai92,P)=aoi92)il + PAyi92)) for cos^^2 e (35°, 80°) , ^ 180°; (5.9) 

where the ± sign in front of PAy{() term is subject to the Madison convention 
quoted in Section 12.21 As the number of registered events is proportional to the 
(t(C, P), from Equations 15.81 and 15.91 we expect to register more elastic scatterings 
during the spin down mode, since A^^ (6^2) is positive for the quoted 92 ranges [28]. 
Starting with the unidentified spin orientations A and B, we divided the available 
range of ^2 into 10 bins, each of 4° width. For each bin and both spin orientations 
the numbers of elastic scatterings (6*2) and (6*2) have been determined and 
exemplary results for the experiment performed at Q = 10 MeV are presented in 
Table O 



92 







(^2) 


Nb (^^2) 


40 


200054 ± 836 


353972 ± 1106 


44 


159864 ± 748 


275061 ± 983 


48 


128808 ± 679 


214802 ± 874 


52 


103183 ± 612 


162223 ± 761 


56 


71068 ± 513 


105549 ± 616 


60 


65849 ± 493 


90496 ± 575 


64 


61430 ± 469 


76328 ± 520 


68 


53476 ± 437 


61810 ± 467 


72 


37984 ± 364 


42024 ± 384 


76 


14784 ± 232 


15694 ± 236 



Table 5.1: Number of elastic scatterings (^2) and (^2) for the unidentified spin 
orientations A and B as a function of tiie CM scattering angle 62 as measured during 
the experiment with pbeam = 2010 MeV/c (Q = 10 MeV). 

Having a look at the numbers in Table 15.11 one can notice that for each bin of 
the scattering angle ^2 we have: 

NBi92)>NAi92), (5.10) 

therefore one can identify mode A with the spin up orientation and mode B 
with spin down. The difference between (6*2) and (6*2) is undoubtfuUy 
statistically significant. The same method was used for identification of the spin 
up and down modes for the experiment at Q = 36 MeV. 
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5.3 Calculation of the relative luminosity 

A method of the determination of the relative luminosity is presented. A detection 
subsystem used for this purpose is described. Systematical errors of this method are 
evaluated. 

In order to determine the relative luminosity Lrei of Equation 12.221 the scat- 
tering in the polarisation plane (see page [6]) had been used, as in this plane the 
differential cross section for any nuclear reaction induced by the strong interaction 
does not depend on the magnitude of beam polarisation. This is a consequence 
of the fact, that the parity in strong interactions is conserved. For detailed expla- 
nation the reader is referred to Appendix O 

Therefore, whenever in the experiment we are restricted to the scattering in 
the polarisation plane the result of the measurement should only depend on the 
polar angle 9, and is independent of the beam polarisation. Thus, the number of 
reactions - n{t) - registered in the polarisation plane within a time interval t may 
be used as the measure of the integrated luminosity over this time interval: 



In order to determine the relative luminosity, the detector system schematically 
presented in Figure 15.91 had been used. It consists of four scintillator detectors: 
two round-shaped, placed vertically in front of the COSY-11 magnet, and two 
square-shaped located horizontally between the magnet coils. The location of 
the scintillators has been chosen in order to conform the kinematic conditions 
of the elastically scattered protons in the polarisation plane. The coincidence 
rate originates mainly from the elastic proton-proton scattering, which constitutes 
about 75% [100] of all reactions. The rest of about 25% comes from multibody 
reactions. The number of coincidences of the following type: 



was measured. 

Denoting by and rii numbers of coincidences integrated over the time of the 
measurement during the cycles with spin up and down, respectively, we get from 
Equation 12.221 and relation 15. lit 




(5.11) 



{sCl^sc2)y {sC:i^sCi). 



(5.12) 



(5.13) 



'rel 



The numbers of coincidences n^ and nj^ and also the relative luminosities Lrei 
for runs at Q = 10 MeV and Q = 36 MeV have been determined using this method 
and are presented in Table 15. 2[ 
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Figure 5.9: Schematic view of the detection system dedicated for the relative luminosity 
determination, as mounted at the COSY-11 section of the COSY ring, (a) front view, 
(b) side view. Beam is circulating along the z-axis and target is along the y-axis. 



Q[MeV] 


n| 




Lrel 


T corrected 
^rel 


10 


5732570 ± 2394 


5942438 ± 2438 


0.96468 ± 0.00056 


0.98468 


36 


5657318 ± 2378 


5874621 ± 2424 


0.96301 ± 0.00057 


0.98301 



Table 5.2: Relative luminosities, and numbers needed for their calculation for runs 
at excess energies Q = 10 and Q = 36 MeV. L'^^^^'^^^'^ are the values of the relative 
luminosities, corrected for the shift of the target, as explained in the text below. The 
errors indicated in the table are the statistical uncertainties only. The systematic errors 
are evaluated in the text. 



The detectors were installed centrally around the nominal target position 
(Ax = 0), however one should keep in mind that the Ax found in Section [5.1.41 
equals 0.7 cm. This introduces false asymmetries into the calculations of Lrei, and 
so, the values of L^ei have to be corrected for these false asymmetries. In order to 
do so, Monte-Carlo simulations have been performed, using the GEANT3-based 
code [101], containing the geometry of the luminosity system from Figure 15.91 
and taking into account the phase-space distribution of the events modified by 
the differential cross sections [102] and analysing powers for the pp — > pp reac- 
tion [28]. In the simulations the real values of the polarisations (quoted in the 
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following section) were taken into account. Subsequently, the same number of 
events have been generated for spin up and down modes, and indeed, it has been 
found that the 7 mm shift introduces false asymmetry that changes the value of 
the L^ei from L^ei = 1 (for the symmetric system) to L^^i = 0.98. Therefore the 
value of AL = 1 — 0.98 = 0.02 has to be added to the Lrei in order to obtain the 
real values of the luminosities, which in Table [5^2] we called Zj'^°Y'^'^^<^d. _ 

It appears from results given in Table 15.21 that the corrected relative lumi- 
nosities differ from 1 by circa 2%, although the spin was flipped from cycle to 
cycle. The overall integrated luminosity during spin down cycles is greater than 
the luminosity for spin up. This apparent inconsistency may be explained by the 
fact that the measurements were always started with the spin down orientation, 
however the breaks in the data taking (due to the target regeneration, beam opti- 
mization and other events aiming in improvement of the quality of measurement) 
happened accidentally with the same probability for both spin adjustment^. 

Statistical uncertainties were calculated according to the rule of the error prop- 
agation applied to Equation 15.131 



where 



-(Lre^) = \l {^) ' + ' a(n,)^ (5.14) 



^(^TU)) = (5-15) 



are the statistical uncertainties of After inserting the partial derivative and 
relations 15.151 into Equation 15.141 the formula simplifies to: 




a{Lrel) = Lrel\ \ • (5.16) 



After these remarks on the statistical uncertainties we will discuss the system- 
atic errors. The main source in the systematic error of the relative luminosity 
determination originates from the uncertainty of the position of the scintillator 
detectors from Figure 15.91 A possible non-symmetrical adjustment of the detec- 
tors with respect to the polarisation plane would lead to the false asymmetries and 
hence to the error of the Lrei determination. In order to estimate the systematic 



^Indeed, during the experiment there were circa rt = 16 breaks in the operation of the 
COSY accelerator. Therefore, in average, there were around 8 cycles with spin down more than 
the ones with spin up. The whole time of measurement for experiment at Q = 10 MeV lasted 
T = 829201 s. The time period of spin up or spin down cycle for the run at Q = 10 MeV was 
equal to t = 1200 s. The rough estimation gives the total number of cycles (spin up and down 
altogether) equal k = T/t « 690 - out of which, we assume, there were 349 cycles with spin 
down and 341 with spin up. Therefore the relative difference in the time of measurement with 
spin down and spin up mode equals to circa 2%. The same inference is valid for the experiment 
at the excess energy Q = 36 MeV. 



5.4 Beam polarisation 



53 



error the conservative assumption has been made that the centres of the scintil- 
lators from Figure 15.91 were shifted from the nominal position by Ar = 4 mm. 
In order to estimate the systematic uncertainty of L^e/, again we have performed 
Monte-Carlo simulations, taking into account all the factors quoted when dis- 
cussing the influence of the Ax on the Lrei- Numerical evaluation brings the 
value of 1% for the systematic uncertainty of Lrei- 



5.4 Beam polarisation 

The methods of the determination of beam polarisation along with the used detection 
systems are described. Ah the methods are basicahy based on the asymmetry measure- 
ments for the pp ppr] reaction. 

During two measurements of the analysing power at different excess energies, 
three independent methods have been used in order to determine the beam polar- 
isation. In the run at the excess energy of Q = 10 MeV the COSY-11 polarimeter 
has been used as the main equipment to extract information about the value of 
the polarisation degree. During this run in parallel the polarimeter of the COSY 
team [103] has also been used as the auxiliary detector to monitor the polarisation 
and to verify the information from the COSY-11 polarimeter. 

During the measurements performed at the excess energy of Q = 36 MeV 
we managed to get access to the EDDA polarimeter [103], commonly used by 
many experimental groups to perform the exact determination of the degree of 
polarisation at the COSY accelerator. 

All three methods of the polarisation monitoring will be described in this 
section. 



5.4.1 Measurement at the excess energy of Q = 10 MeV 

The methods for the determination of the degree of polarisation by means of the internal 
COSY polarimeter as well as the COSY-11 polarimeter will be described. 



The internal COSY beam polarimeter 

One of the polarimeters used for the polarisation monitoring in this experiment 
was a polarimeter of the COSY team [103] . Subsequently we will refer to it as to 
the COSY polarimeter. A schematic view of this setup is presented in Fig. 15.101 
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Figure 5.10: Sciiematic view of the COSY polarimeter. For the short-cuts' explanation 
see text. Figure is adapted from [103]. 



Protons from the polarised beam are scattered on the protons from an internal 
CH2 fiber target and subsequently they reach the detectors presented in Fig. I5.10[ 
There is a kinematical correlation between the angles of the forward (O/) and 
backward (Gf,) scattered protons, namely: 

2,771 

tan e f^iab tan Qb,iab = 777—77 ? > (5-17) 

Tp + 27npC^ 

where m^ denotes the proton mass, while Tp is proton's kinetic energy. The geom- 
etry of this polarimeter was designed such that due to this kinematical restriction 
the polarimeter can operate within the beam energy range from 300 MeV up to 
several GeV. 

The detector setup may be rotated via remote control around the beam axis by 
A(f) = 290^ in order to eliminate the false asymmetries. For the true asymmetry 
measurements, the A0 ranging from 23° (for 0/ = 11°) to 13" (for ©/ = 19'^) 
has been used. For the dimensions and other details of the detector geometry the 
reader is referred to [103]. 

The principle of the determination of the degree of polarisation by means 
of this setup is based on the asymmetry (e) measurements of the quasi-free pp 
elastic scattering realized by using a CH2 fiber target. All the detectors presented 
in Figure IF. 101 are scintillator detectors. The beam line is the symmetry axis of the 
detector setup. FWi, FW2, and their symmetric mirrors constitute the forward 
detector to register the fast proton, while the recoil protons are detected in the 
conjugate backward arm, consisting of four detectors BWi, BW4, and a larger 
common detector BW5 which helps to get rid of the accidental coincidences and 
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give the trigger signal along with the FWi and FW2. Thus, the trigger conditions 
to be fulfilled for calculation of the asymmetry are: 

FWl^FW2^BW^, (5.18) 

for left and right arm of the polarimeter. 

The asymmetry e(0): 

,.0^ _ iv+(e)-iv-(e) 

'^^^ - iv+(e) + iv_(e) ^^-^^^ 

is calculated using the geometrical averages N-t defined as follows: 



iv+ = V^i^i = c'(i + PAy{e^)), (5.20) 

and 



N_ = ^nInI = C(l - PAy{e^)), (5.21) 

where C is the constant depending on the detector efficiencies and the relative 
luminosity during spin up and down modes. The arrows denote the spin mode of 
the beam particles. 

Thus, the degree of polarisation may be calculated as: 

where ^4(0) is the analysing power for the proton-proton elastic scattering. Please 
note, that the degree of polarisation calculated in this way is independent of the 
efficiencies of the detectors. 

In order to correct the above formula for false asymmetries, originating from 
the beam misalignment, the following formula has been applied [104]: 

P ^ ^(e)-e-(9) 

A(0)i-e(e)e'(e)' ^^-^^^ 

where e'(B) stands for the false asymmetries, calculated from equation I5.19[ for 
a measurement with an unpolarised proton beam. The results of polarisation 
measurement with the COSY polarimeter [105] are presented in Figure 15.141 in 
the following section. Here, we would only like to mention that the main source of 
the systematic error in the calculations of P arises from the beam misalignment 
and this has been estimated to be less than 4% [105]. 

It's worth mentioning that the construction of this polarimeter allows to cal- 
culate the polarisations for spin up and down separately. It occurs that during the 
experiment the discrepancy between both polarisations was less than ±3% [105]. 
In order to compare the polarisation values obtained with this equipment with the 
results of the COSY-11 polarimeter for time periods when the COSY polarimeter 
has been in operation the average values of the polarisations during spin up and 
down cycles have been calculated. 
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COSY-11 polarimeter 

The second method devoted to the evaluation of the degree of polarisation made 
use of the COSY-11 detector setup. This method is based on the asymmetry 
measurement for the elastic pp ^ pp process. The COSY-11 detectors that served 
for this purpose are presented in Figure 15.51 of Section 15.1.41 

The elastically scattered events have been identified based on the constraints 
given by the two body kinematics. The relation [5. 171 for 6i and Q2 puts constraints 
on the event distribution, and according to that relation there is a correlation 
between the hit position in SI - x^i - and the index of the silicon pad that gave 
a signal. This correlation may be seen in Figure [5.111 a as an enhancement in the 
density of the event distribution. The events lying outside the correlation curve 
originate from inelastic reactions, which occasionally fulfill the trigger conditions 
for the elastic scattering: 

Tez„, = ^1 A ^4. (5.24) 

The background subtraction is necessary for the evaluation of A^_|_(0) and A^_(0) 
- the numbers of elastically scattered events during cycles with spin up and down, 
respectively. 
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Figure 5.11: (a) Correlation plot for the proton-proton elastic scattering at the beam 
momentum pteam = 2010 MeV/c. For the explanation see text, (b) Exemplary spectrum 
of the distance of the events to the correlation curve of Figure (a) for the arbitrarily 
chosen angle @CM = 47''. The dotted line denotes the approximated background limit. 



To facilitate an easy background subtraction, the correlation curve has been 
divided into three ranges as shown in Figure 15. Ili a and in each range it was 
approximated by a straight line. Event distribution visible in this figure has been 
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projected along the correlation line, and an exemplary spectrum, for Qcm = 47°, 
is depicted in figure 15.111 b. Next, the linear interpolation of the background 
has been performed (see dashed line in Figure IS.lli b). for each range of Qcm 
angle. As it is seen in Figure 15.111 b. the background constitutes only about 5% 
of the signal, and since it is smooth and flat on both sides of the signal peak we 
assumed, conservatively, that the systematic error due to the assumption of the 
linearity of the background is less than 10%. Hence, the overall systematic error 
due to the background subtraction is not greater than 0.5%. Another source of the 
systematic error of the number of events in the individual 0cm ranges originates 
from the possible beam position misalignment, which may have been different for 
spin up and down modes. In our estimations we have assumed conservatively a 
2 mm shift of the beam between spin up and spin down modes, a choice which was 
dictated by the geometrical dimensions of the reaction region [27]. The Monte- 
Carlo simulations performed with this 2 mm shift of the beam revealed the value 
of 0.5% of the systematic error. Therefore, the overall systematic error of the 
number of events in the individual ©cm ranges we estimate to be less than 1%. 

After the background subtraction events lying within 2.5 or less distance from 
the correlation line are regarded to be elastics. Number of elastic events for 
scattering during spin up and down cycles allows to determine the asymmetry 
defined in Equation 15.191 
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Figure 5.12: Relation between the centre-of-mass scattering angle and the position in 
the SI detector for the pp — > pp elastic scattering at pbeam = 2010 MeV/c: (a) Monte- 
Carlo simulations, (b) experiment. Dotted lines are the boundaries of the practical 
range where the polarisation is calculated. The reason for a poor resolution outside this 
region is explained in the text. 



Figure 15.121 depicts the relation between the centre-of-mass scattering angles 
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and the position in the SI detector for the pp pp reaction at the beam mo- 
mentum of Pbeam = 2010 MeV/c. The scattering angles for which both protons 
may still be measured in coincidence by SI and Simon ranges from around 35 to 
circa 75 degrees, as depicted in Figure I5.12[ a. However, the practical partition 
of the CM scattering angles that may be used for the polarisation evaluation is 
squeezed to the range of (37;49) degrees (see Figure IHTT^ b). This is due to the 
strong dependence of the cross section on the scattering CM angle and also due to 
the deterioration of the momentum determination for protons passing the dipole 
only on its edge. The experimental distribution of the Qcm angle as the function 
of the hit position in the SI scintillator is shown in Figure [5.12[ b. 

This range was divided into three sections of 4 degrees width: (37;41), (41;45) 
and (45;49) degrees, choice which was dictated to facilitate an easy comparison 
with the data of the EDDA collaboration. As there were no EDDA measurements 
of the analysing power at the beam momentum Pbeam = 2010 MeV/c, we used 
the average of the pp — > pp analysing powers measured at pteam = 1995 MeV/c 
and Pbeam = 2025 MeV/c [28] (see Figure [5713]) . The obtained average values of 
the analysing power for the pp — pp process at the beam momentum of pteam = 
2010 MeV/c are given in Table 15.31 It is worth noting that values of Ay for 
Pbeam = 2010 McV/c and Pbeam = 2025 MeV/c are about the same within the 
statistical uncertainties. 



i 


0. [°] 


A.y 


1 


(37;41) 


0.385 ± 0.013 


2 


(41;45) 


0.376 ± 0.013 


3 


(45;49) 


0.348 ± 0.013 



Table 5.3: Analysing power values at a beam momentum of Pbeam — 

2010 MeV/c for 

the corresponding ranges of the CM scattering angle. 



For the purpose of the analysis it is enough to know only the value of polar- 
isation averaged over the time of measurement. However, in order to check the 
performance of the used polarimeters we have also performed the differential mea- 
surements of the degree of polarisation. The number of scatterings during spin 
up A+(9) (scattering to the right side with respect to the polarisation plane) and 
spin down A^_(0) (scattering to the left side) has been determined for periods of 
20 cycles, corresponding to circa 1.5 h of measurement. The values of A_|_ and N_ 
have been normalized to the corresponding luminosities. We used the following 
formula to calculate the degree of polarisation: 



P _ Ell ^(0.) 

3 



(5.25) 
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Figure 5.13: Analysing power for the pp — > pp elastic scattering as measured by 
the EDDA collaboration at the beam momenta of Pbeam = 1995 MeV/c and Pbeam = 
2025 MeV/c. In the COSY-11 experiment the beam momentum was set to Pbeam = 
2010 MeV/c, which is the middle of the range of beam momenta presented in this 
figure. 



where 



iv+(e,)-iv_(0.)L, 



rel 



Ay{Qi)N+{Qi) + N,{(S)i)Lrei' 



(5.26) 



The ranges of Bj were chosen as indicated in Table I5.3[ L^e/ are the relative 
luminosities defined in equation I5.14[ calculated separately for each time interval 
of polarisation averaging. The method of L^ei calculation has been described in 
Chapter 15.31 

The confrontation of the polarisation values determined with the COSY and 
COSY-11 polarimeters is shown in Figure 15.141 The open circles are the aver- 
aged values (averaged over spin up and down periods) as obtained by means of 
the COSY polarimeter, whereas the full circles show the results of the COSY-11 
polarimeter. 
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Figure 5.14: The degree of polarisation versus tlie time of the measurement as ob- 
tained during the experiment at the excess energy of Q = 10 MeV. Open circles denote 
results obtained by means of the COSY polarimeter [105], whereas the full points were 
determined using the COSY-11 setup. 

It is worth mentioning that the presentation is only for the sake of the com- 
parison between the results obtained by means of both abovementioned methods. 
Due to the technical problems the COSY polarimeter has been operating only for 
about 1/3 of the period of measurement. Over the most of this time the results 
from both polarimeters are in a good agreement. Low polarisation values in the 
first day of the run was due to the ongoing polarisation development by the ac- 
celerator team. Errors in the Figure 15.141 are the total errors, containing both 
statistical and systematic uncertainties. 

Finally, based on the COSY-11 polarimeter the polarisation value for the full 
time of measurement has been computed and equals: 

P = 0.680 ±0.007, (5.27) 

where the statistical error of the polarisation has been given. The systematic 
error of polarisation determination depends on the uncertainty of the analysing 
power Ay{Qi), the error of the relative luminosity L^ei, and the uncertainty of 
the scattering yields N±{Qi). The systematic uncertainty of the Ay{Qi) equals 
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1.2% [28], the Lrei systematic error has been estimated in Section [5731 to be about 
1%, and the systematic uncertainties of the determination of N±{Qi) is not greater 
than 1%. Applying the error propagation method for the systematic uncertainties 
to Equation 15.261 yields the following formula for the systematic uncertainty of the 
polarisation: 

dP dP dP dP 

= \7^\^Ay + ll^l^L^ei + \j^\AN^ + \^\AN., (5.28) 



with 



dP ^ 1 N+ - LrelN. 

dAy Al^N^ + LreiN.f 
dP _ 1 2N+N_ 

dP 1 2LrelN_ 



dN+ Ay (iV+ + LrelN.) 
dP 1 2LrelN. 



2 ' 



^- (5-29) 

Inserting the partial derivatives 15.291 into Equation 15.281 brings the overall system- 
atic error of polarisation to AP = 8%. 

5.4.2 Measurement at the excess energy of Q = 36 MeV 

The detector setup EDDA and the method of the measurement of the degree of polari- 
sation by means of this aparatus will be described. 

For the polarisation monitoring during this run we made use of the EDDA 
detector setup [28]. This setup has commonly been used at the COSY synchrotron 
as an internal high-energy polarimeter since November 1997{§. 

The schematic view of the EDDA detector is presented in Figure 15.151 Here 
we will only sketch the principle of the polarisation measurement. For more 
information about the detector components and for the explanation of detector's 



^ It is worth mentioning that apart from the monitoring of the polarisation at the COSY 
storage ring, EDDA is a weU operating detector that performed many measurements of the 
PP ~^ PP reactions. The EDDA collaboration measm'ed the excitation functions da / d£l{pp, Qcm) 
for the unpolarised proton-proton elastic scattering [102]. Data has been gathered for 108 
different proton kinetic energies, ranging from 240 MeV up to 2577 MeV, and are available 
at [106]. The EDDA group has also performed the measurements of the excitation function for 
the analysing power A^lpp, Qcm) [28] for 77 different proton beam kinetic energies, ranging from 
436 MeV up to 2492 MeV [107]. The centre-of-mass scattering angles of protons - Qcm - from 
30° up to 90" have been covered. Measurements of the analysing power made use of the polarised 
CH2 fiber target [108]. Recently, the spin correlation coefhcients have been measured [109,110] 
at 34 values of proton's kinetic energy ranging from 504 MeV up to 2493 MeV. 
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operation the reader is referred to [28,111]. The EDDA detector consists of two 
cylindrical double layers that covers Qcm range from 30° to 90*^ for the pp pp 
reaction [28]. The construction of the detector is such that about 82% of the full 
solid angle are covered. The inner part is made of helically wound scintillator 
fibers of 2.5 mm diameter, while the outer part consists of 32 scintillator bars (B) 
surrounded by semirings (R). Target used in the experiment was an atomic beam 
target [108], which can operate in the polarised and unpolarised mode. In the 
measurements of the polarisation for COSY-11 run it has been operating in the 
latter mode. 




mem p 



Figure 5.15: EDDA detector. Figure is adapted from [28]. 



Signals from the outgoing protons are registered in the detectors and sub- 
sequently reconstruction of the reaction vertex is performed with the resolution 
better than 2 mm in all three space directions [28]. Tracks of particles are re- 
constructed and the calculation of kinematical variables is possible. Next, the 
asymmetries e for different 9cm angles can be calculated, using the Formula 15.191 
with corrections for false asymmetries. Polarisation is then given by: 

P = , %.y (5.30) 
Ay cos(0) 

where Ay is the analysing power for the pp ^ pp reaction, determined in previous 
experiments with the EDDA detector setup, and is the azimuthal position of 
the hit in the detector. 

As the acceptance of the EDDA detector for registering events scattered to 
the left and to the right side is large, it was possible to extract the information of 
the polarisations for spin up and spin down (P^ and P^, respectively) separately. 
It was found that [112] 



P^ = 0.642 ±0.004 ±0.008, 
P^ = 0.684 ±0.004 ±0.008, 



(5.31) 
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where the statistical and systematic errors are given, respectively. Thus, the 
average polarisation degree equals: 

P = 0.663 ±0.003 ±0.008. (5.32) 

In our analysis we used this averaged value due to the reasons explained in Sec- 
tion 12. 4[ Determination of the polarisation degree for both spin orientations by 
the EDDA collaboration allowed us to estimate the error we made by this aver- 
aging. Similarly as during the measurement at Q = 10 MeV [105] this amounts 
to about ± 3% only. 

5.5 Determination of the background free pro- 
duction rates for the pp pprj reactions 

We show the way of the determination of the production yields used for reckoning the 
analysing power. Different methods at both excess energies have been used. Differences 
are pointed out in the text. 

In Section Ism we presented the missing mass spectra for the pp —>■ pprj reaction 
at two excess energies: Q = 10 and 36 MeV. The histograms in Fig 15.81 include 
all events measured for both beam spin orientations. In order to determine Ay 
as a function of the centre-of-mass polar angle of the 1] meson emission - - we 
have to divide the events according to the 6*^ range and consider events measured 
during spin up and down modes separately. Optimizing the statistics and the 
expected shape of the analysing power function, the range of the centre-of-mass 
polar angle of the t] meson emission for both excess energies has been divided into 
four bins, and the results were integrated over the remaining four kinematical 
variables introduced in Section [2l3l 

Here we will describe the methods of background subtraction from the missing 
mass spectra determined for both excess energies. We shall start with the descrip- 
tion of the procedure for the excess energy Q = 36 MeV, where the shape of the 
background could be parameterized by a polynomial function of the second or- 
der. Then we will describe the more sophisticated methodology of the background 
subtraction at Q = 10 MeV. 

5.5.1 Q = 36 MeV 

In order to extract the production rates the missing mass spectra have been de- 
termined as a function of the cosine of the centre-of-mass emission angle of the t] 
meson, for spin up and down modes separately. The values of Ay will be derived 
in accordance with the Madison convention quoted in Section 12.21 Subsequently 
we will apply the notation defined in Section 12.41 with N]_ and A^i denoting the 
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production yields to the left side with respect to the polarisation plane during 
spin up and down modes, respectively. 

An exemplary spectrum of the missing mass for spin up mode and 
cos^,, G [0.5, 1] is presented in Figure [5?T6l Next, we made the assumption that 
the backgrouncill shape may be described by a polynomial of second order, and 
the function 

Fbackgr{fnfn) = d + h mm + c mm'^ , (5.33) 

where mm stands for the missing mass, has been fitted to the experimental his- 
tograms in the ranges outside the t] peak with a, b, and c treated as free parame- 
ters. 
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Figure 5.16: A histogram of missing mass for the pp pprj reaction at Q 
Spectrum for spin up mode and cosOfj G [0.5, 1] is shown. 



36 MeV. 



Afterwards, the background evaluated in such a way has been subtracted from 
the experimental spectra for each histogram separately. The resulting background- 
free histograms of missing mass are presented in Figures [5. 17[ a-h as full circles. 

Due to low acceptance of the COSY-11 system for the production of the r] 
meson in the backward directions in the centre-of-mass system at the excess energy 
of Q = 36 MeV (missing mass bin for cos 6*,, e [—1,-0.5]), and as a consequence 
of the problems with the background subtraction for this bin, resulting in the 
huge systematical errors, we decided to omit this bin from further analysis and 
determine the analysing powers at this excess energy for cosO^i larger than -0.5. 



^Originating from the multipionic reactions: pp pp2Tr'^, pp — 
pp —> ppTr^TT^Tr~ , pp ppAir'^, pp pp27r*'7r+7r^, and pp pj<27r+27r" 
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PpilT^, 
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Figure 5.17: Missing mass spectra for the frp ppr] reaction at Q = 36 MeV for 
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Subsequently, the simulations of the shape of the missing mass spectrum for 
different ranges of cos 6rj have been done using a program based on the GEANT3 
code [101]. This program contains the exact geometry of the COSY-11 detector 
system as well as the precise map of the magnetic field of the dipole magnet. 
Also the momentum and spatial beam spreads, multiple scattering, and other 
instrumentation effects have been taken into account. Generated events which 
fulfilled conditions equivalent to the experimental trigger have been analysed in 
the same way as the experimental data and corresponding histograms for the 
pp — > ppr] reaction have been determined resulting in three missing mass func- 
tions gpp^pprjlmm, cos Off), where mm denotes the missing mass. Further, to the 
experimental histograms presented in Figure 15.171 the functions 

qUD {mm, cos 6*^) = 7^^^-' (cos 9^)gpp_^pp^{mm, cos 9^) (5.34) 

have been fitted, for spin up and down separately. Note that the functions 
gpp^ppri{mm, cos 6*^) describing the missing mass shape are the same for spin up 
and down modes, and vary only with the cos 6'^. On the other hand, the free 
parameters of the fit - 7^^-'-''(cos 6',,) - might be different, depending on the spin 
orientation. The number of identified rj mesons for each cos 6r^ bin and for spin up 
and down modes have been calculated from the determined 7^^-''''(cos6'^) parame- 
ters, namely: 

N\_Acos6j^)=^^^^\cos6y^) I gpp^ppri{mm, COS 6r,)d{mm). (5.35) 

Jo 

The statistical errors a (A;^(^^))(cos^^) have been estimated based on the formula: 

rmmmax 

a{N^I^_^{cos6r,)) = a{'y^^^\cos6y)) / gpp^ppr^(mm,cos9^)d{mm), (5.36) 

J 

where a{'y^^^\cos6^)) are the estimates of the 7^^^''(cos^^) parameters uncertain- 
ties (standard deviations) determined by means of the MINUIT minimization 
package [113]. Results are presented in Table 



cos 9n 


Nl{cos9^) 


Ai(cose^) 


[-0.5; 0) 
[0;0.5) 
[0.5; 1] 


103 ± 16 
144 ± 16 
259 ± 24 


100 ± 18 
153 ± 18 
296 ± 28 



Table 5.4: Number of identified rj mesons in the pp pprj reaction at the excess 
energy of Q = 36 MeV as the function of the spin orientation and the cosine of the 
centre-of-mass polar angle of the r] meson emission - 
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5.5.2 Q = 10 MeV 

At this excess energy the upper-energetic tale of the rj peak disperses to the edge 
of the kinematical hmit, as presented in Figure 15.181 This makes the identification 
of the upper part of the multipionic background more inacurate, and hence the 
method of the background subtraction similar to the one presented in l5.5.1l would 
be biased by the larger systematical uncertainties. 
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Figure 5.18: A histogram of missing mass squared for the pp pptj reaction at 
Q = 10 MeV for selected cos 0^ range as measured during spin down mode. 



Therefore, to separate the actual production rates from the background, the 
reactions with multipionic production as well as the events with the t] meson 
production have been simulated in order to reconstruct the shape of the missing 
mass spectra. For simulations, the same program as the one described in Sec- 
tion 15.5.11 has been used, but this time apart from the pp pprj reaction also the 
background reactions like pp pp 2n, Svr, An have been generated. 

In order to perform the credible comparison between experiment and the sim- 
ulations, in the Monte-Carlo calculations the geometry of drift chambers, as well 
as the position and geometrical parameters of the target have been fixed to be the 
same as for the experimental data analysis. 

Let the fpp^pp2n{mm'^,ap, Ap, cos 9r,), fpp-,pp3^{mm'^,ap, Ap, cos 9r,), 
/pp^pp47r(^^^, o'p, Ap, COS Or,) , and fpp^ppr,{mim? , ap, Ap, cos 0^) denote the 
generated background functions for double, triple, fourfold pion production 
and for the 7] production, respectively, where mm? stands for the missing mass 
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squared, whereas cXp and Ap are the momentum spread and the deviation of the 
beam momentum from its nominal value. The functions of the type: 

f^^^^mm^, ap, Ap, cos6'^) = q;^^-'^-'(cos^^^) fpp^pp2-K{jnm^ ^ ap, Ap, cos 6^) + (5.37) 

7^(-'')(cos6'^) /pp_pp4^(mm^ dp, Ap, COS 6*^) + 
5'f^-''^(cos6'^) fpp^pp^{mm^ , dp, Ap, cos 6*^), 

with q;(cos6'^), . . . 5(cos6'^) as free parameters have been fitted to the exper- 
imental spectra of the missing mass squared using the MINUIT [113] mini- 
mization package. For Q = 10 MeV also the bin for cosft^ G [—1,-0.5] 
has been considered in the analysis, as the statistic was sufficient to ex- 
tract the numbers of events also in this particular bin. The fit has been 
performed simultaneously to 8 histograms of missing mass squared, each for 
different spin-cos range configuration. It is important to point out that 
for a given cos 6*^ bin the shapes of the generated missing mass histograms, 
i.e. the functions /pp^pp27r(^"^^, Cp, Ap, cos 6*,,), fpp^pp3T,{mm'^ , ap, Ap, cos^^), 
fpp^ppiTrirnrn? , ap, Ap, cos^^) and fpp^ppr^imrn? , ap, Ap, cos 6'^) depend only on the 
width of the momentum distribution ap and the shift from the optimal beam mo- 
mentum Ap. The of fit has been minimized as a function of 34 parameters: 
amplitudes a{6i),...,5{6i) of the generated background and signal reactions (4 for 
each histogram) and 2 parameters responsible for the spread and the absolute 
value of the beam momentum. However, it is worth noting that altogether the 
histograms contained 480 points. 

In Fig. 15.191 the missing masses for the individual cos On subranges are shown. 
Full dots denote the experimental data, the shaded parts of the histograms depict 
the multipionic background, whereas the solid line represents the best fit of the 
Monte-Carlo data to the experimental spectra. The minimum value of fo^' this 
fit divided by the number of degrees of freedom was determined to be 1.6. 

Numbers of 7] mesons A^^^(_-)(cosft^) for the individual ranges of cos^^ for spin 
up and down cycles have been calculated as: 

N^^_^{cos6rf) = 6^^^\cos6.ri) / fpp^ppnimm^ ,ap, Ap, COS 6n)d{mm^). 

J 

(5.38) 

The statistical errors a (A^^y^^(cosg^)) have b een estimated based on the formula 
analogous to Equation I5.36t 



<^(^+(''^)(cos6'.^)) = a{b^'^^\cosQT^) \ jpp^pprl^mm^ , a p, Ap, cos Qri)d{mm 







(5.39) 

where a{p'^^\cosQr^) are again the estimates of the standard deviations of the 
^T(i)(^(.Qg^^^ parameters. The obtained number of r] mesons per cos 6n range are 
quoted in Table 1531 
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Figure 5.19: Missing mass spectra for different ranges of cos^^ for spin up (a-d) 
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cos 9r, 


Nl{cos9r,) 


iVi (cos e^) 


[-i;-o.5) 

[-0.5; 0) 
[0;0.5) 
[0.5; 1] 


306 ± 27 
267 ± 22 
198 ± 18 
279 ± 23 


250 ± 26 
260 ± 24 
208 ± 19 
286 ± 25 



Table 5.5: Number of registered mesons as the function of the cosine of the centre- 
of-mass r] emission angle (0^) and spin orientation of the beam. 



In the table only the statistical uncertainties are presented. The systematic 
errors for the extraction of number of events have been estimated by comparison 
of the values in Table 15.51 with the number of events determined using the other 
method, namely assuming that the background is linear in the range of missing 
mass peak. In such a way for each cos^^ bin we obtained values which differed 
from the numbers quoted in Table 15.51 by no more than 1.5%. This value we 
assigned to the systematic error of N^'^-'')(cos6'r,) at both excess energies. 



5.6 Results — the analysing power 



The results of the analysing power determination are summarized in the table and in 
the figures. 

Having presented the methodology of Lrei, P and A^^j^_-j(cos0i) determination 
we can now equate the analysing power for both excess energies. For this purpose 
we shall exploit Formula I2.23[ which reads: 

A,(cos^,) = (5.40) 

Relative luminosities can be found in Table 15.21 average polarisations have 
been determined in Section 15.41 and finally Tables 15.41 and 15.51 provide the values 
of background free productions rates for Q = 36 and 10 MeV, respectively. The 
analysing powers for the pp — ^ ppr] reaction which were calculated basing on 
these numbers are shown in 3'''^ column of Table 15.61 Values presented in the 
4*^^ and 5**^ column of this table are the mean values over a cos 6^ bin of the 
theoretical predictions for Ay according to the models with the dominance of the 
pseudoscalar and vector meson exchanges, respectively. In the derivation we took 
advantage of the fact that for an isotropic distribution of the differential cross 
section ^ (9) = const, the average analysing power over an angular range 9 is an 
arithmetical average of the analysing powers for the individual 9n subranges of 9 
- for proof see Appendix [Bl 
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Q[MeV] 




Aj^(cos 9n) 








[-l;-0.5) 


0.163 ± 0.099 ± 0.022 


0.133 


-0.170 


10 


[-0.5; 0) 


0.035 ± 0.091 ± 0.012 


0.067 


-0.092 




[0;0.5) 


-0.021 ± 0.095 ± 0.011 


-0.082 


0.092 




[0.5; 1] 


-0.003 ± 0.088 ± 0.009 


-0.144 


0.170 




[-l;-0.5) 




0.147 


-0.001 


36 


[-0.5; 0) 


0.039 ± 0.179 ± 0.012 


0.046 


0.000 




[0;0.5) 


-0.029 ± 0.122 ± 0.010 


-0.062 


0.000 




[0.5; 1] 


-0.084 ± 0.100 ± 0.011 


-0.154 


0.001 



Table 5.6: Analysing power for the pp — > pprj reaction determined at the excess 
energies Q = 10 and 36 MeV. A^*^"'^ and Ay^"^ are the theoretical predictions according 
to the models with the pseudoscalar and vector meson exchange dominance, described 
in Section [3.11 The first quoted error is statistical, whereas the second one is systematic. 
The estimation of the systematic errors will be given in the next section. 



Statistical uncertainties of Ay were calculated according to the rule of error 
propagation, and the following formula has been applied for their determination: 



a{Ay) = 
with 



\ \dLre 



<J^ (Lrel) + 



dAy 

dP 



(P) + 



dL 



rel 



2N\N 



dAy 

dNl 



T Ari 



P 



Ni + UelN 



2 ' 



dAy 



(iVi) 
(5.41) 



dP 



OA, 



dAy 



1 M - LrelN. 



p2 



2 ' 



2 ' 



2Lre/A^l 



p 



2 • 



(5.42) 



It is worth mentioning, that the partial contributions to the overall statistical error 
from the Lrei and P uncertainties were about two orders of magnitude smaller than 

T I 

the contributions from A^l and uncertainties. 
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Figure 5.20: Analysing power for the pp — > pprj reaction as function of tlie cosine of 
polar emission angle of the r] meson in the centre-of-mass system for Q = 10 MeV (a) 
and Q = 36 MeV (b). Full lines are the predictions based on the pseudoscalar meson 
exchange model [11] whereas the dotted lines represent the results of the calculations 
based on the vector meson exchange [14]. In the right panel the predictions of the vector 
meson exchange dominance model are consistent with zero. Error bars in both panels 
of the figure show the statistical uncertainties only. Systematic errors are estimated in 
next section. 



The analysing power from Table along with statistical errors and the models 
predictions for Q = 10 and 36 MeV are depicted in Figure 15.201 a and I5.20[ b. 
respectively. 



5.7 Systematic uncertainties of the analysing 
power 



We estimate the systematic uncertainties for the analysing power at both excess energies. 



The physical quantities needed to calculate the analysing power Ay, according 
to the Equation I5.40[ are the average beam polarisation P, the relative luminosi- 
ties between spin up and down modes Lreu and the production rates N^'^^^ for 
spin up and down, respectively. 

The systematic uncertainty of polarisation for Q = 10 MeV has been estimated 
in Section [5.4. II and equals AP(Q = 10 MeV) = 0.055. The same quantity for the 
excess energy of Q = 36 MeV was found to be AP(Q = 36 MeV) = 0.008 [112]. 



5.7 Systematic uncertainties of the analysing power 
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Differences in the systematic uncertainties between the measurements are due to 
the fact that different detector setups were used in the determination of the degree 
of polarisation. The systematic uncertainty of the relative luminosity has been 
found in Section [5l3] and equals ALrei = 1% for both excess energies. Finally, the 
systematic error of the determination of the production rates has been estimated 
in Section [5.5.21 and amounts to AN^^^^^-j = 1.5%, the same value for both excess 
energies. 

The values quoted above have been used to determine the overall systematic 
uncertainty of A^,. The contributions from different error sources were added 
analogously to Formula 15.281 exchanging P< — > Ay, N_(_ < — > and N„ < — > Ni. 
This procedure yielded the values of the systematic uncertainties of A^^ quoted in 
Table EH 

At this point it is worth mentioning that the beam misaligment between spin 
up and down modes has no significant influence upon the determination of N_^\\ 
values. Monte-Carlo simulations analogous to the ones presented in Sections 15.31 
and l5. 4.1] have been performed and the systematic error of Nl[^l^ originating from 

the beam misaligment was found to be negligible in comparison with AN_^^_y 

It has also been confirmed that the systematic uncertainty of A^^_(_) originating 
from different binnings of histograms from Figures [5. 1 71 and 15 . 1 91 is negligible with 
comparison to the quoted 1.5% uncertainty originating from the determination of 
A^' using different background models. 



Chapter 6 



Interpretation of the 
experimental results 

6.1 Test of the predictions for A^^ 

Using the statistical inference it is shown that the derived values of the analysing power 
disagree with the predictions of the vector meson exchange models at a significance level 
of a = 0.006. The predictions of the pseudoscalar meson exchange dominance model 
are in line with the data at a significance level of of 0.81. The values of the amplitude 
of the analysing power are extracted at both excess energies. 

In order to determine the statistical significance of our results we have com- 
pared the predictions of the models with the experimentally determined values of 
the analysing power for both excess energies. For both theoretical hypotheses a 
value of the have been calculated according to the formula: 

with a (AyY) denoting the experimental uncertainty of the value and i enu- 
merating the points for the excess energies of Q = 10 and 36 MeV. Altogether 
there are 7 experimental points determined in both experiments. The A^"'^'^'' are 
the values of the analysing power calculated according to the pseudoscalar and 
vector meson exchange models, averaged over the same angular divisions as was 
performed for the experimental data. The value of the analysing power for an 
angular range 6 has been taken as the mean value of the individual analysing 
powers predicted for this range (see Appendix [B]) . These values along with the 
experimental values of Ay^[ are quoted in Table 15. 6[ 

We have obtained the values of equal to 3.78 and 19.32 for the pseudoscalar 
and vector meson exchange models, respectively. Taking into account that there 
are 7 degrees of freedom, the reduced value of the fo^' the pseudoscalar meson 
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exchange model was found to be Xpsd ~ 0.54, which corresponds to a significance 
level apsd = 0.81, whereas for the vector meson exchange model xlec — 2.76 
resulting in a significance level a^ec = 0.006 0. 

Assuming that the predictions of the vector meson exchange model regarding 
the shape of the angular dependence of the analysing power are correct, we have 
also performed the t^st in order to determine the amplitude of the angular 
dependence of the analysing power. As mentioned in Section 13.11 the angular 
distribution of the analysing power may be parametrized with the following equa- 
tion: 

Ay{Q,e,) = A~^(Q)sin20„ (6.2) 

where the amplitude A'^°'^'^'^'^{Q) is a function of the excess energy Q, and its 
energy dependence is shown in Figure 13.71 as the dotted line. A fit of the 
abovementioned function to the experimental results has been performed at 
both excess energies, with A™'^^''"'^'^ as the only free parameter. We have found 
p^max,vec ^ _ q q^i ± Q Q^g fo^ Q = 10 MeV, and A;;^"^'^^'^ = - 0.081 ± 0.091 
for Q = 36 MeV. The best fit functions of Equation 16 21 to the experimental data 
are presented in Figure 16.11 
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Figure 6.1: Fit of the function Ay{Q,9^) = A^"'^'^'^'^ {Q) sin26n to the experimental 
results of analysing power for Q = 10 MeV (a) and Q = 36 MeV (b). The purpose of 
the fit was to extract the amphtude ^"■^''"'^'^ for both excess energies. 

Similar studies have been performed for the pseudoscalar meson exchange 
model. Although authors of this model do not give an explicit analytical pre- 
scription of the analysing power's angular dependence, comparison of the predic- 
tions for different excess energies leads to the conclusion that they do not differ 

^The significance levels as the functions of the values have been taken from [114] 



6.1 Test of the predictions for A 
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much in shape. This may be seen in Figure I6.2[ a. where the theoretical line for 
Q = 10 MeV is compared with the theoretical line for Q = 25 MeV. The latter 
curve has been normalized such that the difference between the predictions for 
these two excess energies is at its smallest. One can see in Figure [6^ a. that the 
differences are indeed much smaller than the experimental accuracy of the deter- 
mined analysing power. Therefore, we may assume that the pseudoscalar meson 
exchange model's predictions with respect to the analysing power have the same 
shape for all excess energies within the close-to-threshold region, and only the 
amplitude of this function varies for different excess energies. Thus, by analogy 
to Equation 16.21 we may write that 

Ay{Q,e,) = A^'^^'P^''{Q)f{e,), (6.3) 

where /(6'^) is the function presented in Figure [6^ a. 

Based on this assumption we have performed a fit of the predicted theoret- 
ical functions for Q = 10 and 36 MeV to the experimental data at correspond- 
ing excess energies, minimising the value with varied as a free pa- 
rameter. The amplitudes A^'^^'P'"^ = - 0.074 ± 0.062 for Q = 10 MeV and 
p^max,psd ^ _ Q Qgg ^ Q^Qg Q = 36 McV havc bccu derived. 




cos9 [ MeV ] cos(e„) cos(9„) 

a) b) c) 



Figure 6.2: (a) Comparison of the shape of the angular distribution of A^^ based on 
the pseudoscalar meson exchange dominance model for arbitrarily chosen Q = 10 MeV 
and 25 MeV. Predictions at Q = 25 MeV were normalized to those of Q = 10 MeV. 
The numerical values of Ay{Q,9>ri) have been made available by the authors of the 
pseudoscalar meson exchange dominance model [11, 115]. (b) Fit of the Ay function 
for Q = 10 MeV to the experimental data, (c) Similar fit, but for excess energy of 
Q = 36 MeV. 

The values of the analysing power amplitudes together with the theoretical 
predictions are depicted in Figure [6731 a for the pseudoscalar meson exchange model 
and in Figure I6.3[ b for the vector meson dominance model. 
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Figure 6.3: Theoretical predictions for the amphtudes of the analysing power angular 
dependence in the close-to-threshold region confronted with the amplitudes determined 
in the experiments at Q = 10 and Q = 36 MeV. The solid line in Figure (a) shows the 
prediction based on the pseudoscalar meson exchange dominance model [11], while the 
dotted line in Figure (b) corresponds to the prediction of the vector meson dominance 
model of reference [14]. 

One can see in this figure that although the amplitude of the vector meson ex- 
change model for the excess energy of Q = 36 MeV lies within about one standard 
deviation from the experimental result, the predicted A^"'^ for Q = 10 MeV is 
4.3 a away from the data point. In the pseudoscalar meson exchange dominance 
model predictions of A^"-^ lie within one standard deviation from the experimental 
data. 

Another observation that can be made from the results shown in Figure 16.31 is 
that the extracted A^"^^ for different models are of about the same value. This is 
due to the fact that the pseudoscalar and vector meson exchange models predic- 
tions do not differ much in the shapes of the angular dependencies of the analysing 
power. 

6.2 Mechanism of the 77 meson production in 
nucleon-nucleon coUisions 

The most probable mechanism of the rj meson production in proton-proton collisions is 
presented. 




As pointed out in Section 13.1.11 it has been established by the theoretical 
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analysis of the data from the total cross section measurements for the pp 
pprj reaction, that the resonant current constitutes the main contribution to the 
production mechanism of the 7] meson in proton-proton collisions in the close- 
to-threshold energy region. It appears that amongst the available resonances 
the Sii(1535) resonance plays the most important role as an intermediate state. 
Excitation of the proton to this resonance may proceed by either one of the light 
pseudoscalar or vector meson exchanges. However, as far as only the data on the 
total cross section are concerned an univocal statement cannot be made on which 
out of the spectrum of possible mesons gives rise to the Sii(1535) resonance. Here, 
the possible particles to be exchanged are the vr, rj, p, and u mesons. 

Some limitations to the models have been made by the measurements of the 
total cross sections for the quasi- free pn — ^ pnr] reaction [13]. It has been found 
that the a{pn — >■ pnri)/a{pp — >■ ppr]) ratio is rather constant in the wide excess 
energy range from 16 to 109 MeV, and equals about 6.5. From such strong isospin 
dependence of the production amplitudes it has been deduced that the excitation 
of the Sii(1535) resonance proceeds via the exchange of the isovector mesons, 
hence this discovery discarded the rj and u mesons as possible intermediate par- 
ticles. Therefore, there remained two possible particles relevant to excite the Sn 
resonance: the vr meson, belonging to the light pseudoscalars meson's nonet and 
the p meson - a member of the light vector meson's nonet. 



N N 




1^1 (1535) 

71 




N N 

Figure 6.4: The most probable mechanism of r] meson production in the nucleon- 
nucleon collisions: excitation of the nucleon to the Sii(1535) resonance via the exchange 
of a TT meson, and its further decay into a nucleon-r/ system. 

We claim to solve this ambiguity by the measurement of the angular depen- 
dence of the analysing power presented in this dissertation. The results are pre- 
sented in Table 1531 of Chapter [521 as well as in Figure 15.201 The analysis of the 
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results allows us to reject the predictions of the vector meson dominance model 
at a significance level of 0.006. Predictions of the pseudoscalar meson exchange 
model are in line with data at the significance level of 0.81. This fact, together with 
the abovementioned presented inference [116] indicate that in nucleon-nucleon 
collisions the t] meson is produced predominantly by the exchange of 
a 71 meson. The most probable process of t] meson production is presented in 
Figure 16. 4[ 

6.3 Ay and the final state of the ppr/ system 

Based on the experimental data we infer on the final state of the rj meson in the ppr] 
system. 

The analysing power determined in this thesis is within the statistical accuracy 
consistent with zero, at both excess energies. This result implies that the rj me- 
son is predominantly produced in the s-wave in the close-to-threshold 
region [59]. This observation is in agreement with the results of the analysing 
power measurements performed by the DISTO collaboration [15] where, interest- 
ingly, in the far-from-threshold energy region Ay was also found to be consistent 
with zero within one standard deviation. 



Chapter 7 
Perspectives 



Some ideas on extending the experiments presented in this thesis are pointed out. An 
experiment for resolving the final state partial waves of the pprj system is proposed. 
It is mentioned that the measurements with the new WASA-at-COSY facility would 
significantly increase the statistics. 

The improvement of the statistics concerning the data on the analysing power 
for the pp — >■ pprj reaction would be possible by means of the measurement of 
this observable with the recently brought to operation An detector WASA-at- 
COSY [69]. Due to the installation of a pellet target, high luminosities for the 
experiments with polarised proton beams are expected to be achieved, and are 
estimated to be at the order of 6-10^^cm~^s~^. This would yield around 20000 rj 
events per day measured at an excess energy of 10 MeV and about 70000 events at 
an excess energy of 36 MeV. Therefore, a one-week measurement of the analysing 
power using the WASA-at-COSY detection setup would yield the result with circa 
20 times better statistical significance than measurements reported in this dis- 
sertation. The letter of intent for such an experiment has already been published 
by the COSY-11 collaboration [117] and is awaiting realisation. 

Another interesting experiment would be the measurement of the spin corre- 
lation coefficients Cxx for the pp pprj reaction. As proposed by Nakayama et 
al. [59], measurements of this observable would be helpful to extract the contri- 
butions of the individual partial waves for the pp pprj reaction in the close- 
to-threshold energy region. Namely, restricting to the rj meson production in 
the s-wav^ and final proton-proton state in the 5* and P waves there are only 
three partial wave^ that can possibly contribute to the pp — > pprj reaction [27]: 

-'^Evaluation valid for the experiment performed at the excess energy of Q = 36 MeV 
^Which is confirmed by the results presented in this thesis. 

•^Here the following notation of the partial waves has been used: ^^'^^Lj — >^'^ '^^ L'j,l, with 
S, L, and J denoting the total spin, orbital angular momentum, and total angular momentum 
of two protons in the initial state. The prime values are corresponding quantities in the final 
state. I stands for the orbital momentum of the rj meson with respect to the pair of protons. 
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^Pq Sqs, ^Sq PqS, and -^^ P2S. Denoting the amplitudes of these 
transitions by a, (3, and 7, respectively, yields [59]: 



da 
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(7.1) 



with k and p standing for the relative momentum in the final and initial proton- 
proton system, respectively, and x = k-p. From the above equation it is clear that 
the quantity ^ {C^x + 1) depends only on the amplitude a. Therefore measuring 
the spin correlation functions in the region where the differential cross sections 
for the pp —>■ ppr] reaction are known would provide the model-independent in- 
formation about the magnitude of the ^Pq Sqs transition. In a similar way, 
the combination ^ {Cxx — 1) different than zero would indicate the presence of 
higher partial waves in the proton-proton final state. Authors of [59] also provide 
the predictions for the Cxx angular dependence, which are presented in Figure EH 
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Figure 7.1: Predictions of the spin correlation coefficient Cxx for the pp — > pprj 
reaction as a function of the final proton angle in the overall centre-of-mass system at 
the excess energy of Q = 41 MeV. The meaning of the curves is explained in the text. 



The dashed line in Figure mi assumes ^Pq -^^SqS transition only, leading to 
the constant value of Cxx equal to 1. The dash-dotted line apart from ^Pq -^^Sqs 



The values of the angular momenta are expressed in the spectroscopic notation: L = S, P, D,. 
and 1 = s, p, d,. . . 
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transition takes into account also the transition ^5*0 ^^PqS. In this case the 
authors also obtained a constant value of C^x ~ —0.3. Inclusion of the —^^P2S 
contribution results in the dotted line. Taking into account the higher partial 
waves would end up with the angular dependence of Cxx being given by the solid 
curve in Figure EH 

The abovementioned experiment has been proposed by the COSY-11 collabo- 
ration [117] to be performed on the new WASA-at-COSY detector setup. However 
there are many technical challenges to be fulfilled before the determination of this 
measurement will become possible. 



Chapter 8 
Summary 



We have presented the theoretical background, the method of measurement and 
the results of experiments aiming in determination of the analysing power for the 
pp — > pprj reaction in the close-to-threshold energy regime. Measurements have 
been performed utilizing the polarised proton beam of the COSY accelerator, the 
cluster jet target dehvering jets of H2 molecules, and the COSY-11 experimental 
facility used to register and identify the reaction products. For the 77 meson iden- 
tification the missing mass method has been applied. The monitoring of relative 
luminosity has been realized with a dedicated detection subsystem, presented in 
Section [5^ measuring the differences in the numbers of reactions taking place in 
the polarisation plane during the spin up and down modes. The degree of polari- 
sation has been determined by means of a series of measurements of asymmetries 
for the pp ^ pp process, utilising three independent polarimeters. Experiments 
have been performed at beam momenta of pbeam = 2.010 and 2.085 GeV/c, which 
for the pp — ^ pprj reaction correspond to the excess energies of Q = 10 and 36 MeV, 
respectively. Results of the data analysis are summarized in Section 15.61 

For the first time ever it was possible to experimentally pin down the domi- 
nating production mechanism of the 7] meson in nucleon-nucleon collisions. Our 
results indicate, at a significance level of 0.81, that the vr meson is an intermediate 
boson exchanged between colliding nucleons in order to excite one of them to the 
resonant state Sii(1535). In the latter part of the process, this baryonic resonance 
deexcites with emission of a nucleon and the rj meson, as presented in Figure 16. 4[ 
It is important to note, that not only the measurements presented in this disserta- 
tion contributed to this finding, but also many hitherto performed investigations 
by various experimental groups have been important in the understanding of this 
process [2-10,12,13]. 

We have also shown that the predictions of the analysing power angular de- 
pendence according to the vector meson exchange model [14], where the p meson 
plays the most important role as an intermediate particle exciting the nucleon 
to the Sii state, disagree with the experimental data at a significance level of 
(^vec = 0.006. This result is intuitively clear and understood, as the branching 
ratio of the Sii(1535) decay into Nvr channel is about 45%, which is an order of 
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magnitude larger than the branching ratio for a decay into a Np pair, measured to 
be at circa 4% [99] . Moreover, the couphng constant for the NNn vertex is larger 
than the corresponding coupling constants for the NNuu, NNp, NNrj vertices. 

One should, however, keep in mind that the interferences in the exchange of 
different types of mesons are not excluded and should be studied theoretically and 
experimentally by the measurement of further spin observables. 

The analysing power of the pp pprj reaction for both excess energies studied 
in this work was found to be consistent with zero within one standard deviation. 
This may suggest that the rj meson is predominantly produced in the s-wave, 
an observation which is in agreement with the results of the analysing power 
measurements performed by the DISTO collaboration [15] where, interestingly, in 
the far-from-threshold energy region the Ay was found to be also consistent with 
zero within one standard deviation. 

The results of this dissertation might be helpful in revisiting the theoretical 
models, in the sense that they may provide a new input with respect to the cou- 
pling constants used for modeling the rj meson production in hadronic interaction, 
and the range of the reaction. 
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Appendix A 

Pseudoscalar and vector mesons 



Definitions of the pseudoscalar and vector mesons as well as differences between isoscalar 
and isovector mesons are given. The structures of the mesons building up the pseu- 
doscalar and vector mesons nonet are given. 

According to QCD, mesons are bound states of quark q and antiquark q' . The 
quarks q and q' may be the same or different. As quarks and antiquarks are the 
spin 1/2 particles, they may form triplet states (|t) with the total intrinsic spin 
J = 1, and singlet states (ti) with J = 0. 

By convention, each quark is assigned positive parity and each antiquark has 
negative parity. If L denotes the orbital angular momentum of the qq' pair, then 
the parity of the meson built out of this pair equals P = (— 1)^"*"^. 

From the three lightest quarks - u, d, and s - whose properties are quoted in 
Table lATTj nine possible qq' combinations can be built. This set of SU(3) mesons 
includes the octet and a singlet states, which can be schematically written as: 

3®3 = 8©1. (A.l) 

Let us assume that qq' are the ground state combinations of quark-antiquark 
pairs with the relative angular momentum L = 0. This implies the parity P of a 
such constructed meson equals —1. With this condition, the mesons with internal 
spin J = are called pseudoscalar mesons and the ones with J = 1 are called vector 



Quark 


B - baryon 


Q - electric 


I3 - isospin 


S - strangeness 




number 


charge 


third component 




u 


1/3 


+2/3 


+1/2 





d 


1/3 


-1/3 


-1/2 





s 


1/3 


-1/3 





-1 



Table A.l: Additive quantum numbers of the SU(3) quarks. By convention, the 
flavours of the quarks (I3 and S) have the same signs as their charges Q. 
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mesons. The ground state pseudoscalar and vector meson nonets are depicted in 
Figure lA.ll 




Figure A.l: The nonet of the ground state (a) pseudoscalar (J^ = 0~) and (b) 
vector (J^ = 1~) mesons. The strangeness S of the meson is plotted versus the third 
component of its isospin I3 . The neutral mesons at the centre of the S-I3 plane are the 
pure mixtures of uu, dd, and ss states. 

States with the same additive quantum numbers, and also the same isospin, 
total internal spin, and parity can mix. This is a consequence of the SU(3) sym- 
metry breaking. The singlet SU(3) state of the pseudoscalar meson nonet - rji - 
corresponds to the following combination of quarks: 

r]i = —={uu + dd + ss) , (A. 2) 

V3 

which mixes with the r^g state, belonging to the SU(3) octet: 

Tjs = —^{uu + dd — 2ss). (A. 3) 

These two states are not the real physical objects. Real mesons which can be 
observed in the experiments - rj and rj' - are the mixtures of these pure SU(3) 
states: 

T] = cos{9psc)r]8 - sin{6psc)r]i, 

V' = sin{epsc)ri& + cos(6'psc)^?i, (A.4) 

where a pseudoscalar mixing angle 9psc = —15.5° has been introduced [118]. 

Similarly, within the vector meson nonet the theoretical SU(3) states ui and 
Us mix with each other. This mesons have the quark content corresponding to 



91 





( iiipitk r*rnn Ml ti pi 1"ioti 

Vc/ LlcXl IV LjVJlliUlllcXljlvJll 


Mas'^ fMpVl 

IVJ-dOO IVIC V 


7r+ 


lid 


139.57 


7r~ 


du 


139.57 


7r° 

K+ 


-k= [uu — dd) 

US 


134.98 
493.68 


K° 


ds 


497.67 


K- 


us 


493.68 




ds 


497.67 


7? 

/ 


A\ (dd + uu) + Bi (ss) 


547.30 


v' 

f 


Ao (dd + uu) + Bo (ss) 


957.78 




OnPiTK pnrn}~tirip^f inn 


Mass fMeVl 


H 


ud 


769.3 


o~ 

r 


du 


769.3 




(^uu — dd) 


769.3 


K*+ 


us 


891.66 




ds 


896.10 


K*~ 


us 


891.66 




ds 


896.10 


UJ 


Ci (^dd + uu) + Di (ss) 


782.57 




D2 [dd + uu) + D2 (ss) 


1019.42 



Table A. 2: Pseudoscalar and vector mesons as the quark-antiquark combinations. 
Masses of the mesons are taken from [120]. 



the rji and r^g states, respectively. The physical states uj and 0, analogous to the 
abovementioned rj and rj' pseudoscalar combinations, are the result of the ui and 
Us mixing with a mixing angle 6yec = 37° [119]. 

The quark structure of the ground state pseudoscalar and vector mesons to- 
gether with their masses are presented in Table IA.2I 



The mesons with the total isospin 1 = 0, like the 1] and rj' mesons within the 
pseudoscalar meson nonet and the and u mesons in the case of the vector meson 
nonet are called the isoscalar mesons. 

By the isovector mesons we refer to either pseudoscalar or vector mesons with 
a total isospin 1 = 1. These are the vr and p mesons, for the pseudoscalar and 
vector meson nonet, respectively. 



Appendix B 

Property of the average analysing 
power 



A proof of the following theorem is presented. 

Theorem 1. For an isotropic distribution of the differential cross section 
^ (9) = const, the average analysing power over an angular range A6 is an 
arithmetical average of the analysing powers for the individual A6i subranges of 



Proof. Without loosing generality, let us consider the production of the rj mesons 
during spin up mode in the plane perpendicular to the polarisation plane. Let 
the Nl (AG) and Nji (A9) denote the number of the rj mesons produced into the 
AG range, symmetrically with respect to the polarisation plane to the left and 
right side, respectively. If P denotes the beam polarisation then, according to 
the Formula l2.19l the averaged beam analysing power for the r] meson production 
into the AO angular range reads: 



AG." 



Ay {A9) 



El, Ay (AO,) 



(B.l) 



n 



Ay (AG) 



1 iVz.(AG)-iVfi(AG) 
PNLiAQ) + NniAQy 



(B.2) 
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Applying the Madison convention defined in Chapter 12.21 and taking into ac- 
count that the production yields A^^^ and are proportional to the cross section 
of Equation 12.131 and under assumption that ^ ( A^) = const we can write: 



Nl (AG) ~ AG (1 + PAy (AO)) , (B.3) 
NniAQ) AQ{l-PAy{AQ)). (B.4) 



Dividing the A0 range into n identical subranges A6i we can rewrite Equa- 
tions EJ and EH 



Nl (AG) ~ J2 (1 + P'^y (^^*))' 

i=l 
n 

Nn (A0) ~ J] Aft, (1 - PAy {AO,)). (B.6) 
1=1 

Putting Equations IB. 51 and IB. 61 into Formula IB. 21 yields: 

, f.r.. ^ 1 ELi Aft. (1 + PAy (Aft)) - ELi (1 - P^y m)) 7^ 

' PEr=iA^^(i+^^y(A^o)+Er=iM(i-PA,(Aft))' ^ 

which after reductions leads to: 



Ay (AWj = =^ — — = — = . (B.8j 

Z^i=iAft r^Aft n 

□ 
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One should mention at this point that the inference presented here does not 
necessarily demands the subranges to be of the same width. This is because 
one can always divide an arbitrarily chosen subrange AOi into m infinitesimally 
small subsubranges, division for which the prove has just been presented. 

For n — > oo the relation B.8 becomes: 

t Ay {9) de 

Ay (Ae) ^ -^^1 ; ^ ^ . (B.9) 



Appendix C 

Parity invariance and the 
analysing power 

It is shown that according to the parity invariance rule the analysing power in the 
polarisation plane equals zero, hence the scattering yields in the polarisation plane may 
be used as an absolute measure of the luminosity. 

The parity transformation is the reflection of the system with respect to the 
origin of the reference frame. Parity is a multiplicative quantum number, which 
is conserved by the strong and electromagnetic interactions. 

Parity transformation may also be represented by a threefold reflection of the 
system in a mirror situated in the y — z, x — z, and x — y planes. This operation 
transforms the right handed (x,y,z) frame into the left handed (x',y',z') frame, as 
depicted in Figure IC.li 



IP 



Figure C.l: Parity transformation. 



Please note that spin, and as a consequence the polarisation vector, is invariant 
under the parity reversal and that these objects transform as the pseudovectors^. 

^T\ie other example of the pseudovector is the angular momentum, which has the same 
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Parity invariance and the analysing power 



Indeed, if we considered the particle with its spin along the y axis we would 
notice that the spin vector changes its direction upon reflecting it in the mirrors 
positioned perpendicularily to x and z axis, but does not change the direction 
when mirrored in the plane perpendicular to the y axis. As a consequence spin 
does not change its direction under the parity transformation. 

Now, without loosing generality, let us consider the scattering of a particle 
with spin along the y axis in the x — z plane, as depicted in Figure IC.2[ a. We 
shall also assume that there is an asymmetry in the left-right scattering in the 
X — z plane, i.e. that the scattering yields A^i and N2 are different. 



Figure C.2: Parity reversal for scattering in the plane perpendicular to the polarisation 
plane. 



The parity reversal transforms the state presented in Figure IC.2[ a into the 
one in Figure IC.2[ b. Please note, that the transformed state of Figure IC.li b 
is physically identical with the initial state from Figure IC.2[ a. as there exist an 
invariant transformation - which is a rotation by angle tt around the y axis - that 
transforms the flnal state into initial one. Naturally, the (x',y',z') frame remains 
a left-handed frame, but the physics of the process does not depend on the choice 
of the reference frame. 

Situation is different when we consider the scattering of a particle with a spin 
along the y axis in the polarisation plane. Again we assume that scattering yields 
A^^i and N2 are different. The initial state and the state resulting from the parity 
transformation are depicted in Figure [UTSl a and IC. 31 b. respectively. 



transformation properties as a spin. 




a) 



b) 
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Figure C.3: Parity reversal for scattering in the polarisation plane. 



In this case no invariant transformation exists that could retransform the final 
state into initial one. We note that the rotation about the y axis does not work in 
this the directions of Ni and N2 after the rotation are reversed. This sit- 

uation is therefore physically not allowed as it would violate the parity invariance 
rule, which holds for strong interactions. 

However, the initial and transformed states would be identical if Ni and N2 
were equal. Therefore, if the studied reaction is invariant under the parity trans- 
formation it cannot result in the asymmetry of yields in the polarisation plane. 
This demands the analysing power Ay in the polarisation plane must be equal to 
zero. 

Therefore the scattering in the polarisation plane does not depend on the 
degree of beam polarisation, and the scattering yields in the polarisation plane 
may be used as an absolute measure of the luminosity. 
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